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ABSTRACT 


This  document  describes  a  theoretical  and  experimental  study  of 
the  interaction  of  high  power,  high  frequency  radio  waves  with  the 
lower  ionosphere.  The  theoretical  calculations  presented  here  show 
that  the  electron  temperature  of  the  ionospheric  plasma  can  be  greatly 
enhanced  when  the  plasma  is  irradiated  by  a  powerful  groundbased  HF 
transmitter  with  an  effective  radiated  power  of  the  order  of  100  MW. 

If  this  plasma  heating  is  maintained  for  times  exceeding  a  few  seconds, 
the  composition  of  the  plasma  can  also  be  altered.  These  temperature 
and  composition  modifications  cause  significant  changes  in  the  plasma 
conductivity  and  wave  absorption  in  the  medium. 

Two  experiments  were  conducted  in  order  to  test  for  the  predicted 
absorption  and  conductivity  modifications:  a  vertical  incidence  pulse 
absorption  experiment  and  a  nonlinear  demodulation  experiment.  Data 
from  the  absorption  experiment  clearly  show  a  large  (9  dB)  increase  in 
wave  absorption  at  2.4  MHz  due  to  a  high  power  (60  MW  ERP)  HF  heating  of 
the  ionosphere.  The  nonlinear  demodulation  experiment  generated  strong 
VLF  radiation  when  the  ionosphere  was  irradiated  by  a  powerful  modulated 
HF  wave.  These  VI^  signals  are  believed  to  be  due  to  HF  heating  induced 
conductivity  modulation  of  the  dynamo  current  system. 


CHAPTER  I 


INTRODUCTION 

1.1  The  Ionospheric  Plasma 

The  ionospheric  plasma  is  a  partially  ionized  region  in  the 
earth's  upper  atmosphere  extending  from  50  km  to  over  1000  km.  It 
consists  of  a  mixture  of  free  electrons,  positive  and  negative  ions 
and  neutral  particles,  and  is  a  transition  layer  between  the  non- 
ionized  lower  atmosphere  and  the  fully  ionized  plasma  of  the  magneto¬ 
sphere.  The  ionosphere  is  formed  and  maintained  by  the  ionization  of 
the  neutral  atmosphere  by  solar  radiations  and  energetic  particles 
such  as  cosmic  rays.  In  the  uppermost  regions  of  the  atmosphere  the 
neutral  particle  density  is  very  small;  hence  there  are  few  particles 
to  ionize.  As  one  moves  to  lower  heights  the  neutral  density  increases 
exponentially,  and  as  a  result  there  is  a  corresponding  increase  in 
the  ionization  density.  However,  the  increase  in  neutral  density  also 
leads  to  increased  attenuation  of  the  ionizing  radiation;  consequently, 
the  ionization  density  eventually  reaches  a  peak,  and  then  decreases 
as  the  ionizing  radiation  is  further  attenuated  with  decreasing  height. 

The  ionosphere  is  usually  divided  into  three  regions  based  on 
the  magnitude  and  structure  of  the  electron  density  distribution. 

These  regions  are  designated  in  order  of  increasing  height  by  the  let¬ 
ters  D  (50  to  90  km) ,  E  (90  to  130  km)  and  F  (>  130  km) ,  respectively. 
The  normal  daytime  electron  density  distribution  has  a  maximum  in  the 
F-region  near  300  km  of  about  10®  cm”3  and  a  secondary  maximum  in  the 
E-region  of  about  105  cm-3.  In  the  D-region  the  electron  density 
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decreases  by  more  than  three  orders  of  magnitude  from  104  cm-3  at 
90  km  to  less  than  101  cm-3  at  50  km.  The  ionospheric  electron  density 
distribution  can  vary  considerably  from  day  to  night  and  during  periods 
of  enhanced  solar  activity. 

In  addition  to  the  ionizing  radiations,  the  physical  properties 
of  the  ionospheric  plasma  such  as  its  temperature,  composition, 
permeability,  and  conductivity  depend  upon  the  motion  of  its  constit¬ 
uent  particles,  which,  in  turn,  depend  upon  the  force  fields  acting 
on  the  plasma.  In  the  ionosphere  these  force  fields  include  gravity, 
the  earth's  magnetic  field,  the  Coulomb  fields  of  the  charged  plasma 
particles,  short  range  force  fields  associated  with  particle  inter¬ 
actions  and  the  electromagnetic  fields  of  radio  waves  which  propagate 
through  the  ionosphere. 

The  focus  of  this  study  will  be  on  the  interaction  between  radio 
wave  fields  and  the  ionospheric  plasma.  However,  in  the  development 
of  this  topic,  it  will  be  necessary  to  consider  the  effects  of  many  of 
the  other  fields  listed  above. 

The  ionosphere  is  an  extremely  interesting  medium  for  the  pur¬ 
pose  of  studying  electromagnetic  wave  propagation.  The  ionosphere,  as 
illustrated  by  its  electron  content,  is  a  nonhomogeneous  medium;  hence 
its  electromagnetic  properties  have  a  spatial  dependence.  Due  to  the 
presence  of  the  earth's  magnetic  field,  the  plasma  exhibits  prefered 
directions  for  wave  propagation  and  is  thus  an  anisotropic  medium. 

The  ionosphere  is  also  both  temporally  and  spatially  dispersive. 
Finally,  and  of  foremost  interest,  nonlinear  electromagnetic  effects 
manifest  themselves  quite  readily  in  the  ionosphere.  An  outstanding 


example  of  this  is  the  well-known  Luxembourg  or  cross  modulation 
effect  (e.g.  Fejer,  1955)  . 

When  an  electromagnetic  wave  propagates  through  a  plasma,  the 
Lorentz  force  due  to  the  wave  fields  alters  the  motion  of  the  elec¬ 
tron  component  of  the  plasma  and,  to  a  lesser  extent,  the  ion  compon¬ 
ent.  If  the  plasma  is  nonrelativistic  (i.e.  one  in  which  the  particle 
velocities  are  much  less  than  the  speed  of  light) ,  the  component  of 
the  Lorentz  force  due  to  the  magnetic  field  of  the  wave  is  much 
smaller  than  that  due  to  the  electric  field.  Hence,  charged  particle 
motion  in  the  plasma  is  perturbed  primarily  by  the  wave's  electric 
field  (Holt  and  Haskell,  1965) .  Now  the  electromagnetic  properties 
of  a  plasma,  such  as  its  permeability,  conductivity,  refractive  index 
and  absorption  coefficient  are  determined  by  plasma  particle  motion 
and  therefore  can  be  altered  by  the  wave  field. 

There  are  several  mechanisms  by  which  the  electric  field  of  a 
radio  wave  can  alter  the  electromagnetic  properties  of  the  ionosphere. 
Two  of  the  most  important  of  these  mechanisms  are  collisional  heating 
of  the  plasma  electrons  (thermal  effect) ,  and  compression  of  the  plasma 
(striction  effect) . 

The  thermal  effect  is  of  principal  importance  in  those  regions  of 
the  ionosphere  where  particle  collisions  are  frequent  (i.e.  the  lower 
ionosphere) .  If  the  electron  mean  free  path  (i.e.  average  distance 
between  collisions)  is  not  too  small,  the  plasma  electrons  can  acquire 
a  significant  amount  of  kinetic  energy  from  the  wave  between  collisions 
as  they  are  accelerated  by  the  electric  field.  When  collisions  between 
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electrons  and  ions  or  neutral  particles  occur,  the  average  energy 
transferred  to  these  heavier  particles  is  small.  The  main  effect  of 
collisions  is  to  randomize  the  directed  electron  motion  imparted  by 
the  wave  field.  As  a  result,  the  average  thermal  energy  of  the 
electrons  (i.e.  electron  temperature)  is  increased  and  the  plasma 
is  said  to  be  "heated."  The  electromagnetic  properties  of  the  plasma 
are  directly  dependent  on  the  electron  energy  and  thus  can  be  modi¬ 
fied  by  radio  wave  heating. 

In  those  regions  of  the  ionosphere  where  particle  collisions 
are  less  frequent  (i.e.  the  F-region) ,  plasma  modifications  due  to 
the  striction  effect  become  important.  In  the  striction  effect,  the 
electric  field  of  a  wave  exerts  pressure  on  the  plasma  electrons 
and  thereby  compresses  the  plasma.  This  results  in  a  local  variation 
in  the  electron  density  and,  since  the  electromagnetic  properties  of 
the  plasma  depend  on  its  electron  content,  there  is  a  corresponding 
variation  in  these  properties  (Gurevich,  1978) . 

In  the  study  of  radio  wave  propagation  in  the  ionosphere,  one 
often  assumes  that  the  electric  field  of  the  wave  is  relatively  weak, 
so  that  the  properties  of  the  plasma  are  not  noticably  affected  by 
the  passage  of  the  wave.  However,  as  the  magnitude  of  the  wave  field 
increases,  the  electron  temperature  and  density  of  the  plasma  can  be 
significantly  modified  by  thermal  and  striction  effects,  and  conse¬ 
quently,  the  electromagnetic  properties  of  the  medium  become  dependent 
on  the  magnitude  of  the  wave  field.  The  propagation  of  radio  waves  in 
the  ionosphere  is,  of  course,  directly  dependent  upon  the  electromag¬ 
netic  properties  of  the  plasma  and  hence  a  strong  wave  can  alter  the 
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propagation  of  all  other  waves  whose  paths  intersect  it.  The  iono¬ 
sphere  thus  supports  the  interaction  of  radio  waves  and  it  is  this 
interaction  which  is  responsible  for  the  nonlinear  phenomena  which 
occur  in  the  ionosphere.  In  particular,  wave  interaction  invalidates 
the  use  of  the  superposition  principle  which  characterizes  linear 
media  and  leads  to  a  nonlinear  relationship  between  the  electric 
field  and  the  electric  displacement  and  current  densities  in  the  plasma. 

1 . 2  Ionospheric  Modification 

Ionospheric  modification,  in  the  broadest  sense,  is  an  attempt 
to  modify  the  physical  properties  of  the  ionospheric  plasma  in  a 
controlled  manner.  This  may  be  accomplished  by  the  in  situ  release  of 
chemically  reactive  gases  or  plasma  clouds,  by  bombarding  the  iono¬ 
sphere  with  particle  beams  from  accelerators  aboard  spacecraft,  by 
VLF  stimulation  of  plasma  instabilities,  or  by  irradiating  the  iono¬ 
sphere  with  high  power  radio  waves  (usually  HF  or  higher  frequency) 
from  ground  based  transmitters  (Walker,  1979) .  This  study  is  concerned 
with  the  last  of  these  techniques,  radio  wave  modification,  and  shall 
focus  on  nonlinear  electromagnetic  effects  arising  in  the  lower 
ionosphere  due  to  plasma  heating  (i.e.  the  thermal  effect) . 

Plasma  modification  experiments  employing  high  power  radio  waves 
have  been  ongoing  since  the  late  sixties  at  Platteville,  Colorado 
(Utlaut,  1970),  at  Arecibo,  Puerto  Rico  (Gordon  et  al.,  1971)  and  at 
several  locations  in  the  U.S.S.R.  (Gurevich  and  Shlyuger,  1975).  The 
motivation  for  these  experiments  was  provided  by  theoretical  studies 
(Ginzburg  and  Gurevich,  1960;  Farley,  1963;  Meltz  and  LeLevier,  1970) 
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which  showed  that  the  deviative  absorption  of  energy  from  HF  waves  , 

propagating  through  the  F-region  could  substantially  modify  the 

local  electron  temperature  and  density.  The  results  of  F-region 

heating  experiments  conducted  during  the  seventies  support  this 

theory  and  have  revealed  many  unanticipated  physical  phenomena 

including  artificial  spread  F,  field  aligned  scattering  and  plasma 

line  enhancements.  The  heated  region  also  acts  as  a  significant 

radio  scatterer  at  frequencies  as  high  as  the  UHF  band  and  is  thus 

potentially  useful  for  telecommunications  purposes  (Utlaut,  1975)  . 

Before  any  of  the  energy  radiated  by  a  ground  based  transmitter 
can  be  deposited  in  the  F-region,  it  must  first  traverse  the  lower 
regions  of  the  ionosphere  where  the  large  neutral  densities  and 
hence  high  electron-neutral  collision  frequency  can  result  in  signif¬ 
icant  nondeviative  absorption.  If  the  transmitted  wave  field  is  large, 
the  absorbed  wave  energy  can  significantly  heat  the  plasma.  Theoret- 

* 

ical  studies  by  Holway  and  Meltz  (1973)  and  Meltz  et  al.  (1974)  based 
upon  current  estimates  of  electron  energy  loss  due  to  collisions  pre¬ 
dict  large  electron  temperature  enhancements  in  the  D  and  E  regions 
when  the  plasma  is  irradiated  by  high  power  radio  waves  at  frequencies 
near  or  above  the  local  electron  gyrofrequency  ( ~  1  MHz) .  For  an 
effective  radiated  power  (ERP)  of  100  MW,  a  factor  of  five  or  more 
increase  in  electron  temperature  is  anticipated.  A  temperature  change 
of  this  magnitude  can  significantly  alter  the  electromagnetic  properties 


of  the  plasma,  resulting  in  a  number  of  interesting  nonlinear  phenomena. 


One  of  the  principal  results  of  plasma  heating  by  a  high  power 
radio  wave,  often  called  a  "heating"  wave,  is  the  modification  of 
wave  absorption.  Wave  absorption  in  a  plasma  is  directly  related  to 
the  frequency  of  electron  collisions  with  other  plasma  particles, 
and  the  frequency  of  these  collisions  increases  as  the  plasma  is 
heated.  A  heating  wave  can  therefore  vary  the  amplitude  of  all  other 
waves  that  intersect  its  path  by  altering  the  absorption  of  these  waves 
In  particular,  if  the  heating  wave  is  amplitude  modulated,  its  modu¬ 
lation  can  be  transferred  to  another  wave  via  the  ionosphere  (cross 
modulation) .  Moreover,  if  the  heating  wave  is  of  sufficient  duration, 
it  can  modify  the  absorption  of  itself  (self  absorption,  e.g.  Gurevich, 
1978)  . 

Another  important  class  of  nonlinear  phenomena  arises  from  heating 
induced  conductivity  modifications.  Nonlinear  wave  interactions  via 
the  plasma  conductivity  can  lead  to  the  creation  of  new  waves  at  fre¬ 
quencies  which  are  a  linear  combination  of  those  of  the  interacting 
waves  (nonlinear  mixing)  or,  if  the  heating  wave  is  amplitude  modulated 
at  some  frequency,  the  plasma  conductivity  may  be  varied  at  this 

frequency  and  lead  to  the  generation  of  a  wave  of  frequency  (non¬ 
linear  demodulation,  e.g.  Ginzburg,  1970) . 

If  the  electron  temperature  in  the  ionospheric  plasma  is  main¬ 
tained  at  an  enhanced  level  by  continuous  wave  (CW)  heating,  the 
composition  of  the  plasma  may  be  modified  due  to  the  electron  tempera¬ 
ture  dependence  of  chemical  processes  involving  electrons,  such  as 
ion-electron  recombination  and  electron  attachment  to  neutral  particles 
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An  increase  in  electron  temperature  is  expected  to  decrease  the  rate 
of  ion-electron  recombination,  leading  to  enhancements  of  both  the 
electron  and  positive  ion  densities  in  the  E  and  upper  D  regions 
(Holway  and  Meltz,  1973;  Meltz  et  al.,1974),  while  heating  the  lower 
D-region  may  increase  the  rate  of  electron  attachment  to  neutral 
particles  resulting  in  a  decrease  in  electron  density  and  an  enhance¬ 
ment  of  the  negative  ion  population  (Tomko  et  al.,1980).  In  a  very 
strong  wave  field,  collisional  ionization  of  neutral  particles  by  the 
field  accelerated  electrons  may  also  modify  the  plasma  composition 
(Gurevich  et  al. ,  1977b).  The  electromagnetic  properties  of  the 
ionospheric  plasma  depend  directly  on  the  electron  content  and  thus 
can  be  modified  by  the  above  mechanisms.  In  particular,  heating 
induced  electron  density  variations  could  modulate  the  ionospheric 
dynamo  current  by  changing  the  plasma  conductivity  (Willis  and  Davis, 
1973)  . 

Experimental  evidence  of  heating  induced  plasma  modifications  in 
the  lower  ionosphere  is  at  present  rather  limited,  due  largely  to  a 
lack  of  facilities  with  adequate  heating  power.  The  Platteville  facil¬ 
ity  has  the  capability  of  radiating  100  MW  ERP  (pulsed  or  CW)  over  3 
to  10  MHz  and  is  ideally  suited  for  ionospheric  modification  studies. 
However,  only  a  small  amount  of  preliminary  D-region  observations  were 
made  at  this  facility  during  the  early  seventies  (Utlaut,  1975)  and 
the  facility  was  not  operational  during  much  of  the  later  seventies. 

At  Arecibo,  Showen  (1972)  has  used  the  430  MHz  incoherent  backscatter 
radar  as  a  diagnostic  in  order  to  observe  E-region  heating  at  40  MHz. 


His  measurements  of  backscatter  power  indicate  a  factor  of  two 
enhancement  of  the  electron  temperature  at  90  km.  More  recently. 

Coco  (1979)  has  used  the  430  MHz  radar  as  both  a  heater  and  a 
diagnostic  in  order  to  observe  temperature  modifications  near  100  km. 
D-region  HF  heating  experiments  at  Arecibo  have  been  attempted  using 
a  log  periodic  antenna  mounted  at  the  focus  of  the  Arecibo  Observa¬ 
tory's  305  m  spherical  dish  and  fed  by  a  100  kw  transmitter  (e.g. 

Sulzer  et  al.,  1976).  However,  the  ERP  for  this  arrangement  (about 
15  MW  at  5  MHz)  is  inadequate  for  producing  large  electron  temperature 
modifications.  New  facilities  for  ionospheric  modification  studies 
have  recently  been  completed  at  Arecibo  and  at  Troms#,  Norway  (e.g. 
Kopka  et  al.,  1976),  and  each  has  the  capability  of  radiating  in 
excess  of  100  MW  ERP  (pulsed  or  CW)  over  the  3  to  12  MHz  range.  These 
new  facilities  together  with  the  reactivated  Platteville  facility 
should  soon  supply  a  wealth  of  data  on  radio  wave  modification  of  the 
lower  ionosphere. 

Experimental  results  from  previous  operations  at  Platteville  have 
been  summarized  by  Utlaut  and  Violette  (1974) .  D-region  data  were 
obtained  from  observations  of  vertical  incidence  pulse  (Al)  absorption 
and  cross  modulation.  The  Al  absorption  measurements  were  made  on  a 
2.667  MHz  diagnostic  wave  and  showed  about  a  6  dB  increase  in  attenua¬ 
tion  due  to  heating  at  100  MW  ERP  and  3  to  6  MHz.  Attenuation  of  the 
diagnostic  increased  promptly  ( <  40  ms)  after  turn  on  of  the  heating 
transmitter,  but  recovered  to  its  unheated  level  very  slowly  ( ~  10  min.) 
following  a  long  period  of  CW  heating  (10  min.) .  The  prompt  increase 
in  attenutation  can  be  attributed  to  a  change  in  collision  frequency 


LO 


as  the  electron  temperature  increases,  while  the  long  term  recovery 
is  most  likely  due  to  an  electron  density  modification.  Observations 
of  amplitude  modulation  impressed  on  signals  passing  through  the 
region  heated  by  the  HF  transmitter  (Jones  et  al.,  1972;  Jones, 

1973)  show  modulation  exceeding  8%  at  2.5  MHz,  12%  at  60  kHz  and 
a  few  percent  at  20  kHz  for  heating  at  5.1  MHz  and  50  MW.  Jones 
(1973)  has  estimated  that  a  40%  increase  in  electron  temperature 
between  70  and  80  km  is  necessary  to  reproduce  the  observed  modula¬ 
tion.  Jones  (1973)  also  observed  nonlinear  mixing  effects  at 
Platteville. 

Utlaut  (1975)  has  reported  that  during  a  solar  flare  sudden 
phase  anomolies  (SPA's)  were  observed  in  the  normal  manner  on  VLF  and 
LF  signals  propagating  through  an  unheated  D-region,  but  the  SPA  was 
absent  on  a  60  kHz  signal  passing  through  a  heated  D-region.  Mitra 
(1975a) has  suggested  that  the  absence  of  the  SPA  at  60  kHz  may  be  due 
to  a  change  in  electron  attachment  induced  by  heating,  with  the 
increase  in  electron  density  in  the  lower  D-region  due  to  the  solar 
flare  being  cancelled  by  the  decrease  in  electron  density  due  to 
increased  attachment.  Full  VLF  calculations  by  Tomko  (1978)  support 
this  theory. 

Much  of  the  theoretical  and  experimental  work  on  radio  wave  modi 
fication  of  the  lower  ionosphere  has  been  pioneered  by  Soviet  scien¬ 
tists.  High  power  self  interaction  and  cross  modulation  effects  have 
been  studied  by  Gurevich  and  Shlyuger  (1975),  Krotova  et  al .  (1977), 
and  Gurevich  et  al.  (1977a) ,  while  observations  of  nonlinear  demodula¬ 
tion  effects  have  been  reported  by  Getmantsev  et  al.  (1974)  ,  Budilin 
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et  al.  (1977)  and  Kapustin  et  al.  (1977)  .  Radio  wave  heating  theory 
has  been  extensively  reviewed  and  developed  in  works  by  Gurevich  (1970) 
and  Ginzburg  (1978) . 

Very  recently,  Turunen  et  al.  (1980)  and  Stubbe  (1980)  have 
observed  nonlinear  demodulation  effects  in  Northern  Scandinavia. 

1 . 3  Objectives 

The  objectives  of  this  study  are  to  develop  a  theoretical  model 
which  describes  the  interaction  of  a  high  power  radio  wave  with  the 
complex,  chemically  dominated  plasma  of  the  lower  ionosphere,  and  to 
ascertain  the  validity  of  this  model  by  experimental  observation  of 
nonlinear  electromagnetic  effects  which  result  from  wave-plasma  inter¬ 
actions. 

The  theoretical  foundation  for  this  study  is  given  in  Chapter  II , 
where  Maxwell ' s  equations  are  combined  with  plasma  kinetic  theory  in 
order  to  obtain  a  self  consistent  set  of  equations  which  describe 
wave-plasma  interactions  in  the  ionosphere.  In  Chapter  III,  a  numer- 
cal  procedure  for  the  simultaneous  solution  of  the  electron  energy 
equation  and  the  heating  wave  energy  flux  equation  is  described  and 
implemented  in  order  to  study  electron  temperature  modifications 
arising  from  the  thermal  effect,  and  self -absorption  effects  due  to  the 
action  of  the  modified  plasma  on  the  heating  wave.  In  these  initial 
calculations,  it  is  assumed  that  the  composition  of  the  ionospheric 
plasma,  and,  in  particular,  the  electron  density,  are  not  altered  by 
the  heating  wave.  The  validity  of  this  assumption  is  examined  in 
Chapter  IV,  and  detailed  calculations  of  ion  chemistry  modifications 
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arising  from  radio  wave  heating  are  presented. 

In  Chapter  V  the  electron  temperature  and  density  models  devel¬ 
oped  in  the  previous  chapters  are  used  to  calculate  the  change  in  A1 
absorption  on  a  diagnostic  wave  due  to  plasma  modifications  induced 
by  a  heating  wave.  The  magnitude  and  time  variation  of  A1  absorption 
predicted  by  the  model  is  compared  with  the  experimental  results  of 
A1  absorption  measurements  made  at  the  HF  heating  facility  of  the 
Arecibo  Observatory  by  the  author  and  his  coworkers  (Tomko  et  al., 
1981)  ,  as  well  as  with  the  Platteville  results  which  were  described 
earlier. 

Chapter  VI  presents  the  basic  theory  of  the  generation  of  very 
low  frequency  radiation  by  nonlinear  demodulation  of  high  power  HF 
waves  in  the  ionosphere.  The  electron  temperature  model  developed  in 
Chapter  III  is  used  to  calculate  the  change  in  plasma  conductivity 
due  to  an  amplitude  modulated  HF  wave  which,  in  turn,  is  used  to 
estimate  the  strength  of  ELF  and  VLF  fields  radiated  by  the  ionosphere 
due  to  the  modulation  of  ionospheric  currents.  The  predicted  charac¬ 
teristics  of  the  ELF /VLF  radiation  are  compared  with  the  results  of 
nonlinear  demodulation  experiments  conducted  at  polar  latitudes  (e.g. 
Stubbe,  1980)  and  with  data  from  a  recent  experiment  at  the  Arecibo 
HF  facility. 

The  conclusions  of  this  study  are  summarized  in  Chapter  VII . 

The  research  study  outlined  above  should  prove  to  be  useful  in  a 
number  of  applications.  With  the  increased  availability  of  HF  heating 
facilities,  a  rapid  increase  in  active  experimentation  with  the  iono¬ 
spheric  plasma  can  be  anticipated.  The  ionospheric  model  developed  in 
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this  study  should  prove  to  be  a  valuable  tool  in  the  design  and 
evaluation  of  these  future  HF  heating  experiments.  In  turn,  future 
observational  results  will  undoubtably  lead  to  improvement  of  the 
model  and  a  better  understanding  of  ionospheric  processes.  For 
example,  the  intensity  of  ionospheric  radiation  generated  by  non¬ 
linear  demodulation  of  HF  waves  is  directly  dependent  upon  the  magni¬ 
tude  of  ionospheric  currents  and  their  associated  electric  fields. 

These  radiations  could  thus  act  as  a  highly  sensitive  indicator  of 
natural  changes  in  the  global  current  system.  VLF  and  ELF  radiation 
generated  by  nonlinear  demodulation  may  also  have  potential  applica¬ 
tions  to  submarine  communications. 

The  results  of  this  research  study  will  also  be  useful  in  asses¬ 
sing  the  impact  of  solar  power  satellites  (SPS)  on  the  lower  ionosphere, 
and  on  radio  wave  propagation  in  this  medium.  Proposed  SPS  systems 
(Glaser,  1977)  would  convert  solar  radiations  into  microwave  energy 
which  would  be  transmitted  to  earth  and  converted  into  electrical 
power.  The  SPS  microwave  beam  would  have  a  power  density  as  large  as 
200  W/m2  and  could  thus  cause  significant  electron  heating  in  spite  of 
the  relatively  low  ionospheric  absorption  at  microwave  frequencies. 

The  present  work  extends  the  results  of  previous  SPS  impact  studies 
(e.g.  Perkins  and  Roble,  1978;  Duncan  and  Zinn,  1978;  Holway  et  al., 
1978)  by  taking  into  account  the  complex  chemistry  of  the  lower  iono¬ 
sphere  and  its  dependence  on  electron  temperature. 


CHAPTER  II 


V  X  H(r,  t)  =  D(r,  t)  +  J^r,  t) 

■±  -*■  -►  4-  -► 

V  X  E (r ,  t)  =  -  B(r,  t) 

(2.1) 

(2.2) 

±  ^  ^  >  “► 

(2.3) 

V  •  D(r,  t)  *  P(r,  t) 

(2.4) 

V  •  B (r ,  t)  =  0 

where 

E  -  electric  field  strength 

H  -  magnetic  field  strength, 

-► 

D  -  displacement  density, 

B  -  magnetic  flux  density, 

J  -  current  density,  and 
p'  -  charge  density. 

The  dot  above  a  given  quantity  denotes  partial  differentiation 
with  respect  to  time  (i.e.  9/8t) . 

Strictly  speaking,  the  field  quantities  which  appear  in  Maxwell's 
equations  are  average  values  taken  in  the  statistical  sense  (i.e. 
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averaged  over  an  ensemble  of  equivalent  systems) .  However,  in  plasma 
applications,  one  is  concerned  with  fields  whose  wavelength  is  large 
compared  with  the  mean  distance  between  plasma  particles.  Statistical 
averaging  is  thus  equivalent  to  averaging  over  a  physically  infinit¬ 
esimal  volume,  dr  =  dxdydz,  of  dimensions  which  are  much  larger  than 
the  mean  distance  between  plasma  particles  (i.e.  a  volume  containing 
a  large  number  of  particles) .  Hence  E(r,  t)  is  to  be  interpreted  as 
the  average  value  of  the  electric  microfields  in  the  volume  element 
dr  centered  at  position  r  at  time  t.  The  other  field  quantities  are 

similarly  defined  (Ginzburg,  1970) . 

-►  > 

The  current  density,  J  ,  includes  externally  applied  currents, 

J  ,  such  as  those  associated  with  transmitting  antennae  which  may  be 
EXT 

immersed  in  the  plasma,  and  plasma  currents,  J,  arising  from  the  net 
motion  of  charged  plasma  particles.  Likewise,  the  charge  density, 
p' ,  consists  of  an  externally  applied  part,  p  ,  and  the  plasma  charge 

EXT 

density,  p.  Hence  one  writes 


p' 

(J'  *>  ■  °EXT<J’ 

t) 

f  -+ 

+  p(r, 

t) 

(2.5) 

j' 

<*'  11  ■  3EXT(?- 

t) 

+  J(r, 

t) 

(2.6) 

The  external  charge  and  current  densities  may  generally  be  regarded  as 

■+■ 

given  quantities,  independent  of  the  field  vectors  E,  H,  D  and  B.  The 
charge  and  current  densities  are  related  by  the  charge  continuity 
equation 

$  •  J'it.  t)  =  -  p'(r,  t)  (2.7) 

Maxwell ' s  equations  combined  with  the  charge  continuity  equation 
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amount  to  two  vector  relations  (2.1-2)  between  the  five  vector  unknowns 
E,  H,  D,  B  and  J,  along  with  three  boundary  conditions  (2.3-4  and  2.7)  . 
To  complete  this  system  of  equations,  one  requires  a  set  of  three  vector 
equations  relating  the  displacement,  magnetic  flux  and  current  densities 
in  the  medium  to  the  electric  and  magnetic  field  strengths.  The 
required  set  of  relations  are  called  material  or  constitutive  relations . 

At  an  elementary  level  all  matter  may  be  considered  to  consist  of 
charged  particles;  some  free  (ions,  electrons)  and  others  bound  (neutral 
atoms  or  molecules) .  The  formulation  of  the  constitutive  relations  for 
a  given  material  thus  amounts  to  a  description  of  the  effects  of  the 
motion  of  these  charged  particles  on  the  electromagnetic  field  vectors. 
The  description  of  these  effects  may  be  approached  from  two  different 
viewpoints;  one  macroscopic  and  the  other  microscopic. 

At  the  macroscopic  level  the  effects  on  the  fields  due  to  charge 
motion  are  taken  into  account  by  means  of  a  set  of  constitutive  para¬ 
meters.  The  motion  of  free  charges  is  described  by  means  of  the  con¬ 
ductivity,  a,  the  translational  motion  of  bound  charges  is  described  by 
means  of  the  permittivity,  e,  and  the  rotational  motion  of  charges  is 
described  by  means  of  the  magnetic  permeability,  The  relationship 

between  the  field  vectors  in  the  medium  are  specified  by  the  equations 


D(r,  t)  =  £  E(r ,  t) 

B(r,  t)  =  u  H(r,  t) 

M 

J(r,  t)  =  a  E(r,  t) 

These  simple  relations  are  used  extensively 


(2.8) 

(2.9) 

(2.10) 

in  the  study  of  linear 
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electrodynamics  (e.g.  Jordan  and  Balmain,  1968)  and  apply  to  homo¬ 
geneous,  isotropic  media.  However,  equations  (2.8-10)  are  not,  in 
general,  applicable  in  ionospheric  studies  since  the  plasma  can 
exhibit  nonlinear,  nonhomogeneous  and  anisotropic  characteristics. 

At  the  microscopic  level,  the  electromagnetic  properties  of  a 
medium  are  described  in  terms  of  the  dynamical  behavior  of  the  elemen¬ 
tary  charges  which  comprise  the  medium.  In  this  viewpoint  one  dis¬ 
penses  with  the  concept  of  a  material  medium  per  se ,  and  instead 
considers  the  ensemble  of  plasma  particles  situated  in  free  space. 

In  any  infinitesimal  volume  dr,  it  is  assumed  that  the  number  of 

positive  and  negative  charges  are  equal  so  that  there  is  no  net  charge 
■+  ->■ 

in  dr  (i.e.  p(r,  t)  =0).  The  physical  reasoning  behind  this  assump- 
tion  is  that  if  each  dr  had  a  net  charge,  there  would  be  a  potential 
difference  between  volume  elements  due  to  Coulomb  forces  and  the 
plasma  particles  would  quickly  move  between  volume  elements  in  such  a 
way  as  to  reduce  the  potential  difference  and  restore  neutrality 
(Yeh  and  Liu,  1972) .  In  the  presence  of  an  electromagnetic  field,  the 
charges  in  the  plasma  (free  or  bound)  are  displaced  by  the  Lorentz 
force  due  to  the  field 

F  =  q  (E  +  v  x  B)  (2.11) 

s  s  s 

-*•  — ► 

where  q  is  the  charge  of  the  s-th  type  of  plasma  particle  and  v  (r,  t) 
s  s 

•> 

is  the  velocity  of  some  s-type  particle  in  dr.  This  action  gives  rise 
to  an  average  convection  current 

=  l  %  ns  (  ^s) 
s 


J 


(2.12) 


where  n  (r,  t)  is  the  number  density  of  s-type  particles  in  dr  and 
s 

v  (r,  t)  is  the  average  velocity  of  the  s-type  particles  in  this 
volume.  The  plasma  particles  are  situated  in  free  space;  hence,  the 
permittivity  and  magnetic  permeability  are  those  of  a  vacuum:  £  =  £q, 
y„  =  y  .  One  therefore  has 

M  O 


D  =  £  E  (2.13) 

o 

and 

B  =  p  H.  (2.14) 

o 

Equations  (2.12-14)  are  the  constitutive  relations  in  the  microscopic 
viewpoint . 

Using  the  microscopic  constitutive  relations  one  can  rewrite 
Maxwell ' s  equations  in  the  form 


^  x  h  =  jio£  E  +  y  q  n  (v  \ 
o  L  s  s  \  s/ 
s 

^  x  E  =  -jtou  H 
o 

^  •  D  =  0 
^  •  H  =  0 


(2.15) 

(2.16) 

(2.17) 

(2.18) 


where  it  has  been  assumed  that  the  field  vectors  E  and  H  vary  sinus¬ 
oidally  as  e'lt0t  with  angular  frequency  u)  (i.e.  monochromatic  fields)  , 
and  the  medium  is  assumed  to  be  free  of  externally  applied  sources 

(i.e.  J,,—  =  0,  p__,  =  0)  .  In  order  to  complete  the  system  of  equations 
EXT  EXT 

(2.15-18)  one  must  supplement  them  with  a  set  of  dynamical  equations 


of  motion  for  the  plasma  particles.  The  formulation  of  these  dynamic 
equations  is  described  in  the  following  section. 
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2.2  The  Kinetic  Theory 

The  kinetic  theory  attempts  to  explain  the  physical  properties 

of  a  plasma  by  considering  the  forces  of  interaction  between  the 

component  particles  of  the  plasma  in  the  presence  of  externally  applied 

force  fields.  In  this  theory  the  state  of  the  plasma  is  described  by 

distribution  functions  for  each  type  of  constituent.  The  distribution 

function  f  {r,  v  .  t)  for  the  s-th  particle  type  is  defined  such  that 
s  s 

the  number  of  s-type  particles  in  the  differential  volume  dr  located 
at  position  r  at  time  t  with  velocity  components  in  the  range  to 

-V  -*■  -V  ->■  “V 

v  +  cv  is  given  by  f  (r,  v  ,  t)  dr  dv  where  dv  =  dv  dv  dv  is 
s  s’  s  s  s  s  sx  sy  sz 

a  differential  volume  in  velocity  space  centered  at  Vg,  and  vgx»  VB  ' 

-V 

v  denote  the  cartesian  components  of  v^ .  The  number  density  of 
s-type  particles  at  r  and  t  is  thus  given  by 


uu 

n  (r,  t)  =  (  f  (r,  v  ,  t)  dv 

s  /  s  s  s 


(2.19) 


The  fundamental  equation  governing  the  particle  distribution  func¬ 
tions  is  a  plasma  is-  the  Boltzmann  equation  which  is  of  the  form  (Holt 
and  Haskell,  1965) 


~  f  —  — v  •  V  f  -  R  •  V  f  + 
3t  s  s  s  s  s 


(2.20) 


where  ^  is  the  differential  velocity  operator  defined  by  the  rela¬ 


ys  9  ^  9  .  ~  9  ,  ~ 

v  f  *  -r—  f  x  +  f  y  +  t—  f  z  . 

vs  9v  s  9v  s  9v  s 
sx  sy  sz 


(2.21) 


The  term  R  is  the  force  per  unit  mass  on  s-type  particles  due  to 
s 
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externally  applied  fields.  For  electromagnetic  fields  Rg  is  given 
by  the  Lorentz  force  per  unit  mass 


R  =  q  (E  +  v  x  B  )/m  (2.22) 

s  s  T  s  T  s 

-► 

where  qg  and  m^  are  the  charge  and  mass  of  s-type  particles,  and 
and  Bt  are  the  total  electric  field  and  magnetic  flux  density  acting 

-V  — ► 

on  s-type  particles.  In  this  study  E^  =  E ,  the  electric  field  of  the 
heating  wave  and  Bt  =  y^H  +  Be  where  H  is  the  magnetic  field  of  the 
heating  wave  and  B  is  the  magnetic  flux  density  of  the  earth’s  mag- 
netic  field.  Boltzmann's  equation  states  the  local  variation  in  the 
distribution  function,  is  equal  to  the  variations  in  the 

distribution  function  due  to  streaming  of  particles  in  or  out  of  a 

“f  •+• 

given  volume,  v  •  Vf  ,  externally  applied  forces,  R  •  V  f  ,  and  col- 

S  S  S  Vs 

lisions  of  s-type  particles  with  one  another  as  well  as  with  all 

other  types  of  plasma  particles,  (9f  /9t) 

s  coll 

Particle  collisions  in  a  plasma  are  of  two  basic  types:  long 
range  Coulomb  interactions  between  charged  particles  (e.g.  electron- 


electron,  electron- ion  or  ion-ion  collisions)  and  short  range  inter¬ 
actions  involving  at  least  one  neutral  particle  (e.g.  electron-neutral, 
ion-neutral  or  neutral-neutral  encounters) .  The  former  type  of  col¬ 
lision  is  a  multiple  interaction  in  which  a  charged  particle  continuously 
interacts  with  the  Coulomb  fields  of  a  large  number  of  other  charged 
particles.  In  contrast,  the  latter  type  of  interaction  is  a  localized 
encounter  between  two  particles  (i.e.  a  binary  collision)  .  Based  on 
the  relative  importance  of  each  of  these  types  of  collisions  in  deter¬ 
mining  the  physical  properties  of  the  plasma,  one  can  distinguish 


between  two  different  types  of  plasma:  a  strongly  ionized  plasma  in 


which  Coulomb  collisions  are  dominant,  and  a  weakly  ionized  plasma 
in  which  binary  collisions  overshadow  Coulomb  interactions.  In  the 
lower  ionosphere  the  electron  and  ion  densities  are  much  smaller  than 
the  neutral  density,  and  the  properties  of  the  plasma  are  governed 
primarily  by  electron-neutral  and  ion-neutral  collisions;  hence,  only 
binary  collisions  are  considered  here. 

The  time  rate  of  change  of  the  distribution  function,  f  ,  due  to 

s 

binary  collisions  between  s-type  particles  and  all  other  particles 
can  be  written  as  an  integral  of  the  form  (e.g.  McDaniel,  1964) 


(fs  fk  ~  fs  V  9  qsk(g'  6)  ^  dVk 


(2.23) 


where 

qsk^g'  "  differential  scattering  cross  section  for 

collisions  between  s  and  k  type  particles, 

g  =  lvs  “  vjcl  ~  relative  speed  between  colliding  particles, 

dfi  =  sinQdP  d<{)  -  differential  solid  angle,  and 

0  -  angle  between  v  -  v,  and  v  -  v,  , 
s  k  s  k 

where  the  tilde  above  a  given  quantity  indicates  its  value  following 
the  collision. 

The  collision  integral  includes  all  types  of  binary  processes 
which  occur  in  a  plasma.  Depending  on  the  types  of  particles  involved 

i 

in  the  collision,  these  processes  may  include  elastic  scattering, 
inelastic  encounters  resulting  in  the  excitation  of  rotational,  vibra¬ 
tional  or  electronic  energy  levels  of  one  or  both  of  the  colliding 
partners,  collisional  ionization,  ion-electron  recombination  processes, 


electron  attachment  to  neutral  particles,  charge  transfer  and  many 


other  types  of  chemical  reactions . 

Given  the  distribution  functions  for  a  plasma,  the  mean  value 

of  any  plasma  property,  $  (r,  v  ,  t) ,  associated  with  s-type  particles 

s  s 

(e.g.  velocity,  kinetic  energy,  etc.)  is  given  by 


0  (r,  t) 
s 


1 

n  (r,  t) 
s 


1 


$  (r,  v  ,  t)  f  (r,  v  ,  t)  dv 
s  s  s  s  s 


(2.24) 


—00 

In  particular  the  constitutive  relation  between  the  current  density 
and  the  electromagnetic  field  in  a  plasma  is  given  by 


OO 


J(r,  t)  =  y  q  f  v  f  (r,  v  ,  t,  E,  H)  dv  (2.25) 

L  s  J  s  s  s  s 

s  — oo 

where  the  dependence  of  the  individual  distribution  functions  on  the 
heating  wave  field  vectors,  E  and  H,  is  specified  by  Boltzmann's 
equation. 


2.3  Solutions  of  the  Boltzmann  Equation 

For  a  homogeneous  plasma  on  which  no  external  fields  are  active, 
the  steady  state  solution  of  the  Boltzmann  equation  is  the  well-known 
Maxwellian  distribution  function 


m 


3/2 


s,0 


2TT  K  T 
B  S 


-m_V  •  V  /2  K  T 
,  s  s  s  B 


(2.26) 


where 


T  -  characteristic  temperature  of  the  s-type  particle 
distribution, 


K  =  8.617  x  10-5  eV/  K  -  Boltzmann's  constant,  and 
B 

■+  ^  /  *^  \  , 

V  -  v  -(  v  )  -  thermal  speed  of  s-type  particles. 

S  S  '  s/ 
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if  the  mean  thermal  speed  of  s-type  particles  is  much  larger  than 

2 

their  average  velocity,  one  can  write  V  •  V  s  v  so  that  the 

s  s  s 

Maxwellian  distribution  is  symmetrical  (i.e.  depends  only  on  the 
magnitude  of  v  ) .  This  condition  is  usually  satisfied  for  electrons 
and  ions  in  the  lower  ionosphere  (Gurevich,  1978)  . 

If  the  distribution  function  of  each  particle  type  in  the  plasma 
is  Maxwellian  and  all  of  the  particle  types  are  characterized  by  the 
same  temperature,  then  the  plasma  is  said  to  be  in  a  state  of  thermal 
equilibrium.  In  this  study,  it  will  be  assumed  that,  in  the  absence 
of  heating  wave  fields,  the  plasma  is  in  a  state  of  thermal  equili¬ 
brium  characterized  by  the  neutral  gas  temperature,  T^.  This  is  a 
reasonable  assumption  in  the  lower  ionosphere  (i.e.  <  100  km)  where 
good  thermal  contact  between  the  plasma  particles  is  achieved  due  to 
the  high  density  of  the  neutral  particles  and  the  corresponding  high 
frequency  of  particle  collisions  (Banks  and  Kockarts,  1973)  .  How¬ 
ever,  in  the  upper  ionosphere,  the  neutral  particle  densities  are 
smaller,  particle  collisions  are  less  frequent,  and  the  electron, 
ion  and  neutral  temperatures  can  differ  significantly  (Schunk  and 
Nagy,  1978) . 

For  the  general  case  of  a  plasma  containing  spatial  gradients 
and  acted  upon  by  external  fields,  no  exact  solution  of  the  Boltzmann 
equation  is  known.  However,  one  can  construct  an  approximate  solu¬ 
tion  to  this  equation  by  expanding  the  distribution  function  as  a 
series  of  spherical  harmonics  in  velocity  space  and  them  making  cer¬ 
tain  assumptions  about  the  state  of  the  plasma  in  order  to  make  the 
problem  mathematically  more  tractable. 
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f  , 
s 


The  spherical  harmonic  expansion  of  the  distribution  function, 
in  velocity  space  is  of  the  form  (Brandstatter ,  1963) 

oo  n 


f  (r,  v  ,  t) 
s  s 


l 

n=0 


l 


r  m  / 

f  (r, 
s  ,n 


v  , 
s 


t)  Ym(a, 
n 


6) 


(2.27) 


-*■  in 

where  v  ,  a,  8  are  the  spherical  polar  coordinates  of  v  and  the  Y 
s  s  n 

are  normalized  spherical  harmonic  functions  defined  by 


Ym(a,  6)  =  (-l)m 
n 


(2n+l)  (n-m) ! 


4tt  (n+m)  ! 


*5  .  m  m 

sin  a  d 


P  (cosa)  e 


jmB 


d  (cosa) m 


(2.28) 


where  P  (cosa)  is  a  Legendre  polynomial  of  order  n;  and 
n 


(Ym)*  =  (-1) m  Y~m 
n  n 


(2.29) 


where  (Ym)*  is  the  complex  conjugate  of  Ym.  The  spherical  harmonic 
n  n 

functions  satisfy  the  orthogonality  condition  (Arfken,  1970) 


2TT  tt 

/K1 

0  0 


m 

(a,  8)  Y  2 (a,  6)  sina  da  dB  =  6 
n2 


nl'  n2  mi'  m2 


(2.30) 


where  6 


nl' 


is  the  Kronecker  delta  function, 


1, 

0, 


n^=n 

n^n 


2 

2 


(2.31) 


The  general  procedure  for  solving  Boltzmann's  equation  using  the 

spherical  harmonic  expansion  is  to  substitute  (2.27)  into  the  Boltzmann 

m  * 

equation,  multiply  the  resultant  equation  by  Y^*  and  integrate  over 
sinadadB.  Using  the  orthogonality  relation  for  the  spherical  har¬ 
monics  to  eliminate  the  double  summation  over  n  and  m,  one  generates 


I 


25 

a  series  of  equations  involving  the  functions,  f  m  .  If  one  can 

s  /  n 

solve  these  equations  for  the  f  m  up  to  some  order  n  =  N,  then 

s  ,n 

(2.27)  can  be  used  to  write  an  N-th  order  approximation  to  the 
distribution  function  (e.g.  Gurevich,  1978).  A  modified  form  of 
this  procedure  will  be  used  in  the  following  section  to  study  the 
electron  distribution  function. 


2.4  The  Electron  Distribution  Function 

This  section  outlines  a  first  order  solution  of  the  Boltzmann 

equation  for  the  electron  component  of  a  weakly  ionized  plasma  in 

the  presence  of  a  heating  wave  field  and  a  constant  geomagnetic  field 

of  flux  density,  B  .  The  objective  of  this  analysis  is  to  develop 
E 

— ► 

an  expression  for  the  electron  current  density,  J^,  induced  by  the 
electromagnetic  field  of  the  heating  wave.  Since  one  is  concerned 
only  with  the  field  induced  current,  the  effects  of  gravity  and 
spatial  gradients  will  be  neglected.  It  is  also  assumed  that  the 
plasma  is  nonrelativistic  so  that  the  force  exerted  on  the  plasma  by 
the  magnetic  field  of  the  heating  wave  is  much  smaller  than  that  due 
to  the  corresponding  electric  field,  and  can  thus  be  neglected. 
Finally,  the  electric  field  of  the  heating,  E,  is  assumed  to  vary  har¬ 
monically  as  e^Ut.  Under  these  conditions  the  Boltzmann  equation 
becomes 


9  f  -  e 
9t  e  m 


(E 


+  v  x  B) 
e  E 


V  f 
v  e 


(2.32) 


where  -e  =  q  . 

e 


X 


26 


From  equations  (2.27-28) ,  one  has  to  the  first  order  in  n 


■U) 


h  .  0  (  3  \h  .  ^  jB  1 

f  n  ■  ST  sinot  e  f 
e,0  8ir  e 


*  [*) 


co  sot  f 


+  i,sir)  sina  fe!i 


or  more  compactly 


e ,  1 


(2.33) 


f  =  f  „  +  a  •  f  , 
e  e,0  e  e,l 


(2.34) 


where 


fe,0  ■(«)  C 

^e,l  =  (iiT  )  (fe,l  +  fe,l  ®2  +  fe,l 


S3 


and 


a  =  v  /v 
e  e  e 


j  3  “"i  6 

-sina  eJ  e^  +  cosa  e 2  +  sina  eJ  e^ 


/r 


sr 


A  A 


with  e^,  e3  being  a  set  of  complex  orthonormal  vectors. 

Substituting  (2.34)  into  (2.32)  and  integrating  over  sinadadB, 
one  obtains 


—  f  -  —  t  1  9  ,  2  ■* 

9t  e,0  n>  3v  ^  9v  Ve  e,l 

e  e  e 


where 


2tt  tt 


J.  f 

9t  e,0 I 


)  =  —  f 

1  3t  e,0 I  (2.35) 


coll 


)c.U  '  ^  /i^fe)0011Si"“d“d6 


(2.36) 


0  0 


Similarly,  by  inserting  (2.34)  into  (2.32)  and  integrating  over 
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is  related  to  f  The  corresponding  collision  terms  Of  _/9t) 

e,0  e,l  coll 

and  Of  „/9t)  describe  the  transfer  of  momentum  and  energy, 

e,0  coll 

respectively,  from  electrons  to  the  neutral  gas. 

In  order  to  evaluate  the  collision  terms  (2.36,  2.38)  ,  one  must 

first  evaluate  the  collision  integral  Of  /9t)  ,,  which  is  given 

e  coll 

by  (2.23) .  In  the  weakly  ionized  plasma  of  the  lower  ionosphere, 
electrons  collide  primarily  with  neutral  particles,  principally  N2  and 
O^.  Since  these  particles  are  much  more  massive  than  an  electron, 
one  could  assume  as  a  lowest  order  approximation  that  the  neutral 
particles  are  infinitely  massive  (i.e.  they  are  initially  at  rest 
and  do  not  recoil  as  a  result  of  collisions  with  electrons) .  This 
assumption  implies  that  v  =  v  =  g,  f  =  f  ,  and  since  f  and 

€  ©  1C  K  6fU 

f  depend  only  on  the  magnitude  of  v  ,  one  has  f  =  f  , 

e , 1  G  G|U6#U 

■> 

f  ,  =  f  . .  From  these  conditions  and  the  collision  geometry  given 
e ,  x  e ,  1 

in  Figure  2.1,  one  finds  that 


f  ?  -  f  f  =  (a  -  a  )  •  f  f 

e  k  e  k  e  e  e,l  k 


(2.43) 


=  [sin9  cos<{>  sini|>  +  (l-cos0)]  cos4>  f  f 

6  r  1  1C 


Substituting  (2.43)  into  (2.23)  and  integrating  over  d<J>  and  dv^,  one 
finds  that  the  collision  integral  can  be  written  as 


-if 

3t  e 


coll 


= -I  v  a  •  f  , 
L  ek  e  e ,  1 


(2.44) 


where 


V  ,  (v  )  -2 -n 
ek  e 


7T 

nk  Ve  j 


(1-COS0)  q  ,  (v  ,9)  sin0  d0 
ek  e 


(2.45) 


is  the  electron-neutral  collision  frequency  for  momentum  transfer 
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This  equation  may  be  put  in  a  more  compact  from  by  defining  the 
following  quantities: 


Y  =  -  e  B  /m  oj  (2.51) 

E  e 

Z  =  V  M  (2.52) 

en 

and  U  =  1  -  j  Z  (2.53) 

Rewriting  (2.49)  in  terms  of  these  quantities,  one  has 


-y  -*■  0-^3 

Uf  -  j  Y  X  f  =  -  j  —  E  “  f  „ 
e,l  e,l  m  to  3v  e,0 

e  e 


(2.54) 


Let  x  ,  x2  and  x^  be  a  set  of  orthonormal  unit  vectors  chosen  such 

A  A 

that  Y  is  parallel  to  x^ :  Y  =  Y  x^ .  Writing  the  left  hand  side  of 
(2.54)  in  matrix  one  has 


-y  a  -*■  9 

A  f  -  -  j  E  ~  f  n 

—  e,l  mai  9v  e,0 

e  e 


(2.55) 


where 


A  = 


U  jY  0 

-jY  U  0 

0  0  U 


in  the  above  coordinate  system. 

Solving  (2.55)  by  matrix  inversion,  one  has 


(2.56) 
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(2.57) 


where 


A-1 _ 1 

-  U(U2  -  Y^) 
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(2.58) 
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One  can  now  write  the  electron  current  density  (2.43)  as 

00 

2 


4ir  .  e2  f  -1  3  3  „ 

3  3  m  oj  J  —  Ve3v  e,( 


dv  E 


(2.59) 


From  the  above  equation,  it  is  apparent  that  one  can  write  J  as  a 

e 

conduction  current 


J  =  a  E 
e  — e 


(2.60) 


where  O  is  the  electron  conductivity  matrix  defined  by 


a  = 
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(2.61) 


and 
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(2.63) 


a  „  il  •  _ef.  f  1  v3  3 

L  3  3  m  a)  J  U  Ve3v  e,0 


dv 


(2.64) 


where  0  ,  a  ,  and  O  are  the  transverse.  Hall  and  longitudinal  con- 

T  n  L 


ductivities,  respectively. 


2 . 5  Energy  Transfer 

In  order  to  specify  the  electron  conductivity  completely,  one 
must  solve  equation  (2.35)  for  f  .  As  noted  earlier,  the  collision 

©  f  0 

term  in  (2.35) ,  (3f  _/3t)  , . ,  is  zero  for  the  case  of  electron 

e,0  coll 
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collisions  with  infinitely  massive  neutral  particles?  hence,  there 

is  no  energy  transfer  from  the  electrons  to  the  neutral  particles. 

In  order  to  allow  for  energy  transfer,  the  assumptions  used  in 

deriving  the  collision  integral  (2.44)  shall  be  relaxed  by  letting 

f  (v  )  =*  f  „{v  +  Av  ) .  Since  the  electron  energy  is  governed  by 

e,0  e  e,0  e  e 

f  it  is  again  assumed  that  f  ,  =  f  ,  .  The  assumptions  that 
e,0  e,l  e,l 

the  neutral  particles  are  initially  stationary,  g  =  v  ,  and  that  the 

neutral  particle  distribution  functions,  f^,  are  not  significantly 

altered  by  electron  neutral  collisions,  f,  =  f.  ,  are  also  retained. 

k  k 

Applying  this  set  of  assumptions  to  equation  (2.23)  and  integrating 

over  dv,  ,  one  has 
k 


coJ1  ■  1  \J  I'«,0(V4V  -  fe,0(VK  ®  <2-65> 

k 


-  J 

v  a. 

-*■ 

•  f 

L 

ek  k 

e,l 

k 

Substituting  this 

result 

into  (2.36)  and  integrating  < 

one  finds  that 

/  3  \ 

[  8t  e,o]  ® 

l  n  | 

f  jf  (v  +Av  )  -  f  (v  )  v 

V  I 

'coll 

L  k 

k  J 

i  L  e,D  e  e  e,0  e  J  i 

q  ,  df)  (2.66) 
e  ek 


For  elastic  electron-neutral  collisions,  the  collision  integral 
(2.66)  can  be  written  as  (e.g.  Holt  and  Haskell,  1965) 

I  Jr  f  o)  *  ~~2  V2  (I  G®J  V  )  fv  f  +  (K  T  /m  )  £  l\ 

^3t  eO/  2v2  3veA  e  ek  ek  L  e  eO  B  n  e  3v  eQJ  f 

'  el  Ve  6  t  k  e  J 

(2.67) 
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electrons . 


In  the  lower  ionosphere  electrons  lose  energy  by  inelastic 
collisions  resulting  in  the  excitation  of  rotational  and  vibrational 


energy  levels  in  N2  and  .  For  diatomic  molecules  such  as  N2  and 
02 ,  Gerjouy  and  Stein  (1955)  have  shown  that  the  total  cross  section, 
for  the  excitation  of  rotational  energy  levels  are  given  by 


6. 


J+2 


(2.69) 


where 


=  a  *  «  6  .  „ 

ek  —J  iJ  —  j  ,  J~2 


k  ( J4-2)  (J+l) 

0  (2J+3) (2J+1) 


(4J+6) 

u 


Jn 's 


(2.70) 


k  _  k  J  (J-l) 

-J  “  0  (2J-1) ( 2J+1) 


1  .  "o1 

u 


“o  ■  8,1  2k  *o/15' 


-9 

=  5.29  x  10  cm  -  Bohr  radius, 


-  quadrapole  moment  for  the  k-th  type  diatomic 
molecule  and 

BQ  -  the  corresponding  rotational  constant. 

The  energy  of  the  rotational  quanta  are  given  by 


k 

u.  . 

13 


(4J+6)  B*  6.,  6. 

0  U  3 , 


J+2 


k 

u. 


(4J-2)  B^  6 

0  iJ 


J-2 
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Using  the  cross  sections  (2.69-70) ,  Gurevich  (1978)  has  shown  that 
(2.67)  reduces  to 


I51'60)  tmi?  1[ 


\ot> 


f  n+(KoT  /m  ) 
eO  B  n  e  9v 


(2.71) 


v  v 

provided  that  u»  B  and  K  T  »  B 
0  B  n  o 


k  k 

where  R  k.  (v  )  =  8  B0  ao  nk 

rot  e  - 

m  v 
e  e 


and 


n  =  In. 

k  L  ki 


-2 

In  the  lower  ionosphere,  u~K  T  :2  x  10  eV  while  B  2  =  2.48 

Bn  o 

-4 

x  10  eV  (Herzburg,  1950);  hence,  the  conditions  on  (2.71)  are 
satisfied. 

If  one  writes  (3f  n/dt)  , ,  as  the  sum  of  (2.67)  and  (2.71)  and 
eu  coll 

substitutes  this  sum  into  (2.35),  it  is  found  that  f  is  governed 
by  an  equation  of  the  form 


_3_  =  _1 _ 3_ 

3t  eO  2  9v 
2v  e 
e 


2  e  ■*  ■% 

—  —  E  •  f 

3  m  el 

e 


+  I  G  ,  v  ,  v  f„+KB  n  9 

ek  ek  e  eO  -  eO 

,  >  m  dv 

k  e  e 


(2.72) 


where  G  (v  )  =  G6,1  +  R  k  (v  )/V  (v  ) 
ek  e  ek  rot  e  ek  e 

For  steady  state  conditions  in  which  the  harmonic  perturbations 
to  f  have  been  time  averaged  over  the  period  27T/cu ,  equation  (2.72) 
can  be  reduced  to  (Holt  and  Haskell,  1965) 


i 
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where 


2 

3 


21T 


V 

ek  ek 
[Re  {1} 


0 


V 

e 


K  T 

f  +  -^-2- 
eO  m 


Re  { f  .  }]  dt 
el 


=  0 
(2.73) 


and  Re  { E }  denotes  the  real  part  of  E.  From  (2.57)  one  has 


->•  +  )  f 

f  .  =  r  r-  eO  (2.74) 

el  dv 

e 

±  — -i  0  —1  -> 

where  T  -  -  A  E, 

m  co  — 
e 

therefore 


E  •  f  =  E  •  f  eO  (2.75) 

el  dv 

e 

Substituting  (2.75)  into  (2.73)  and  integrating  over  dv0 ,  one  has 


eO 


m  v 
e  e 


[K  T-  +  2c  E  •  f/3  (  y  G  V  )] 
1  B  n  1  ek  ek' 


(2.76) 


Integrating  once  more  over  dv  ,  one  obtains 

e 

m  v  dv 

_ e  e  e _ 

[  K  T  +  2e  E  •  r/3(  l  G  ,  V  ,  )] 

1  B  n  L  ek  ek  J 

k 

where  the  constant  of  integration,  C,  is  determined  from  the  normal¬ 
ization  condition  (2.19)  : 

00  » 


eO 


C  exp  (  - 


n 

e 


f  f  dv 
l  e  e 

oo 


4tt  f  v2  f  dv 
l  e  eO  e 


(2.78) 
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From  (2.77) ,  it  is  apparent  that  the  symmetrical  part  of  the  distri¬ 


bution  function  is  Maxwellian  only  when  2e  E  •  T/3  (  1  G  .  V  )«  K  T  . 

k  ek  ek  B  n 

In  general,  f ^  is  non-Maxwellian  due  to  the  dependence  of  ?  and 

I  G  V  The  above  analysis  has  included  only  elastic  and  rotational 
k  ek  ek 

energy  transfer  processes;  however,  if  one  includes  vibrational  cooling, 
a  similar  result  is  obtained  (Gurevich,  1978) . 

As  a  lowest  order  approximation,  it  is  assumed  that  the  denomin¬ 
ator  of  the  integrand  in  (2.77)  can  be  written  as  K  T  where  T  is  an 

Be  e 

effective  electron  temperature.  The  symmetric  part  of  the  distribution 
function  is  therefore  approximated  as  a  Maxwellian  distribution 


f  ^  =  n 


/  m  \ 

3/2  j 

1  e  - 

exp  / 

2ir  L  T 

V  B  e/ 

1 

-  m  v 
e  e 


(2.79) 


STe 

and  the  effect  of  the  heating  wave  field  on  the  distribution  function 

is  accounted  for  by  describing  the  variation  in  T^  due  to  the  heating 

wave  field.  The  equation  governing  the  variation  of  Tg  can  be  derived 

2 

from  (2.35)  by  multiplying  this  equation  by  m  v  /2  and  integrating 

e  e 

2 

over  4tt  v  dv  The  resultant  equation  is 
e  e 


where 


3t 


X  n  K  T 
2  e  B  e 


coll 

'e 


coll 

'e 


(2.80) 


(2.81) 


is  the  rate  of  electron  energy  loss  due  to  collisions.  As  in  (2.76), 

the  harmonic  variations  in  T  have  been  time  averaged  over  2tt/u), 

e 
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equation  (2.84),  this  nonlinearity  is  described  in  terms  of  T^  which, 
of  course,  depends  directly  on  the  field  strength  via  (2.80) . 


2 . 6  The  Dispersion  Equation 

In  terms  of  the  kinetic  theory,  the  total  current  density 
induced  by  the  electric  field  of  a  heating  wave  can  be  written  as 


J  (E) 


(E) 


=  l 


v  f  (v  ,  E)  dv 
s  s  s  s 


(2.85) 


where  it  is  assumed  that  f  is  symmetrical  when  E  =  0.  In  the  pre- 

s 

— ►  -> 

vious  two  sections  an  expression  for  J  (E)  has  been  developed.  In 

e 

order  to  completely  specify  J  (E) ,  one  must  also  include  any  field 
induced  ion  currents.  In  order  to  specify  these  currents,  one  could 
proceed  as  in  Section  2.3  to  expand  each  of  the  ion  distribution 
functions  as  a  series  of  spherical  harmonics: 


f.  =  f.„  +  a.  •  f.  +  .  .  .  (2.86) 

i  lO  i  1,1 


The  first  order  solution  of  the  Boltzmann  ecruation  for  the  ionic  com¬ 
ponents  of  a  weakly  ionized  plasma  follows  that  given  in  Section  2.4 
for  the  electrons,  the  primary  difference  being  that,  unlike  the  elec¬ 
tron,  the  mass  of  a  typical  ion  is  of  the  same  order  as  that  of  the 
neutral  particles.  From  elementary  considerations,  Gurevich  (1978)  has 
shown  that  the  ion  currents  are  small  in  comparison  with  the  electron 
current  provided 


(2.87) 
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where  f  =  co/2it 

f  =  - electron  gyrofrequency 

G  2TT  m 

e 

m^  -  mass  of  i-th  ion  type,  and 

V.  -  ion-neutral  momentum  transfer  collision  frequency 
in 

In  the  lower  ionosphere  one  has  m  /m.  <  10“®,  f  =  0.85  -  1.7  MHz 

e  i  G 

(Yeh  and  Lui,  1972)  and  v.  ~  V  (Banks  and  Kockarts,  1973), 

in  en 

hence,  (2.87)  is  satisfied  for  HF  fields.  One  can  therefore  write 

J  (E)  =  J  (E)  . 
e 

Using  this  result  one  can  write  the  first  of  Maxwell's  equation 
(2.15)  as 

^xH  =  D  +  J=  ju>e„  E  +  a  E  (2.88) 

0  — e 

where  D  =  joj£Q  E  -  free  space  displacement  current,  and 

J  =  a  E  -  electron  conduction  current. 

— e 

In  ionospheric  studies,  one  often  equates  the  right  hand  side  of 
(2.88)  to  an  effective  displacement  current 

d'  =  jU)£0  (1  -  j  C^/EgU))  E  (2.89) 

so  that  the  plasma  is  characterized  as  a  dielectric  with  relative  per 
mittivity  tensor 

* 

£  =  (I  -  j  a  /£  <i>)  (2.90) 

~R  — e  0 

In  this  viewpoint  the  currents  induced  in  the  plasma  by  the  electric 
wave  field  are  polarization  currents 

p'  =  jw  P '  =  a  E  (2.91) 

—e 

•+  j, 

where  P  is  the  effective  plasma  polarization  vector  defined  by  the 


relation 
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D  =  tQ  E  +  P' 


(2.92) 


The  above  viewpoint  shall  be  adopted  in  the  following  discussion  of 

the  wave  and  dispersion  equations. 

By  susbtituting  (2.84)  into  (2.90) ,  the  relative  permittivity 

tensor  may  be  written  out  explicity  as 

00 

2 


V“'  V  =  i  -  r 

0  B  e 


A  1  (v  , 
—  e 


uj)  f  _  (v  ,T  )  v  dv  (2.93) 
eO  e  e  e  e 


The  dependence  of  Er  on  T^  leads  to  a  nonlinear  relationship  between 

D  and  E  via  (2.80) ,  while  the  dependence  of  £  on  u)  corresponds  to 

R 

frequency  dispersion.  Frequency  dispersion  arises  from  the  fact  that 

the  effective  polarization,  P' ,  of  the  plasma  in  the  presence  of  a 

harmonic  electromagnetic  field  cannot  respond  instantaneously  to 

changes  in  the  field  due  to  the  inertia  of  the  plasma  charges  and  the 

collisional  forces  to  which  they  are  subject  (Ginzburg,  1970) . 

Following  the  usual  procedure,  the  wave  equation  for  the  electric 
■> 

field  vector,  E,  can  be  derived  by  taking  the  curl  of  (2.16)  and  then 

-V 

applying  (2.88)  to  eliminate  H.  One  obtains 


r,2-*-  2  ,->■ 

V  V  •  E  -  V  E  =  K  £  E 
0  — R 


(2.94) 


U) 


where  K„  =  —  -  free  space  propagation  constant,  and 

0  c 


c  = 


0  0 


speed  of  light  in  vacuo. 


If  one  assumes  that  E  has  the  form  of  a  monochromatic  plane  wave  with 
wave  vector  K,  . 


E  -  Eo  e 


j  (cot  -  K  •  r) 


(2.95) 
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(2.96) 


(2.97) 


(2.98) 


Equation  (2.96)  has  a  nontrivial  solution  if  and  only  if  the  deter¬ 
minant  of  L  vanishes: 

det  (L)  =  0  (2.99) 

Equation  (2.99)  is  the  dispersion  equation  for  a  homogeneous,  aniso¬ 
tropic  plasma.  The  dispersion  equation  describes  the  relationship 
between  the  wave  vector  K  and  the  wave  frequency  in  the  presence  of 
the  plasma.  Since  (w)  is  known,  equation  (2.99)  specified  those 

K 

— ► 

values  of  K  for  which  waves  of  the  form  (2.95)  can  propagate  in  the 
plasma. 

In  order  to  solve  the  dispersion  equation,  one  must  first  eval¬ 
uate  the  integral  in  (2.93)  .  From  equations  (2.52-53)  and  (2.58) , 
one  notes  that  the  velocity  dependence  of  A  ^  arises  from  the  term  U 
which  is  proportional  to  the  electron-neutral  collision  frequency, 

V  .  Thus,  in  order  to  evaluate  (2.93)  one  must  specify  the  velocity 

dependence  of  v  .  In  the  lower  ionosphere  the  neutral  gas  is  about 
en 

78%  and  21%  0^;  hence,  one  can  write  from  (2.50) 

V  (u)  *  n  .  (-2—  Pu*5  [0. 78  a  „  (u)  +  0.21  0  n  (u)l 
en  axrl  m  /  eC>2  1 

'  e 


(2.100) 
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where 


u  =  m  v  /2, 
e  e 


a  (u)  =  2ir 
ek 


(1  -  cos0)  q  ,  (v  ,9)  sin0d0 
ek  e 


-  momentum  transfer  cross  section  for  electron  col¬ 


lisions  with  the  k-th  netural  particle  type,  and 


n  .  -  number  density  of  air. 

air 


The  momentum  transfer  cross  sections  for  N2  and  0 are  given  by 
(Banks  and  Kockarts,  1973) 


O  ,  (u)  =  (18.3  -  7.3  u)  u5  X  10~16  cm2  (2.101) 

eN2 

and 

O  (u)  =  (2.2  +  5.1  u*3)  x  10'16  cm2  (2.102) 

e02 

where  u  is  in  electron  volts.  Figure  2.2  shows  a  plot  of  V  as  a 

en 

function  of  u  based  on  the  above  cross  sections.  For  u  <  1  eV,  one 
has  to  a  good  approximation  that  v  is  proportional  to  u.  Budden 
(1965)  has  evaluated  (2.93)  for  this  case  and  has  shown  that 


(G1  +  £2)/2 

j(ex  -  e2)/2 

0 

-j(e!  "  £2}/2 

(e1  +  e2)/2 

0 

where 

ei  =  1  “  Z  [a0  C3/2  (a0}  +  j  2  C5/2  (<v] 
in  w 

£2  =  1  -  I  [ax  C3/2  <ax>  +  j  I  C5/2  (ax}] 

e .  *  l  -  §•  f  z-1c  ,_(Z-1)  +  j  f  c  (z_1)l 

3  ZmL  m  3/2  m  2  5/2  m  J 


(2.103) 


ELECTRON  ENERGY  (eV) 

Electron-neutral  collision  frequency  as  a  function  of  electron  energy. 
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Cl.  =  (1  +  Y)/Z 
0  in 


a  =  (1  -  Y)/z 

x  m 


2  ,  2 

X  =  u>  M 
e 


/  2  \  *5 

/n  e  \ 

03  =  [  — -  ]  -  angular  electron  plasma  frequency 

e  I  m  e  ' 


e  0 


Z  =  V  /03 
m  m 


v  =  v  (K  T  )  -  electron-neutral  momentum  transfer 
m  en  B  e 


collision  frequency  evaluated  at  the 
most  probably  electron  energy 


and  (Dingle  et  al. ,  1957;  Hara,  1963) 


If  Z* 

C  (X)  =  —7-  — 5s - r 

p  p:  J  (E2  *  x2) 


e  ^  d£ . 


Equation  (2.103)  applies  to  the  coordinate  system  illustrated  in 
Figure  2.3,  which  was  used  in  defining  A  (see  2.56) .  If  one  chooses 
the  wave  vector  to  be  in  the  x^,  x^  plane,  then  one  can  easily  show 


that 


r  K2cos26 


K  •  K  E  -  K  •  EK  = 


0  -K  cos©  sin0 
K2  0 


-K2cos0sin0  0  K2sin20 


(2.104) 


where  E  =  E1  x  +  E2  x2  +  E3  x3  and  0  is  the  angle  between  Y  and  K. 

By  substituting  (2.103-104)  into  (2.99)  and  solving  for  K  (e.g. 
Yeh  and  Liu,  1972) ,  one  finds  that 


K  =  K0  n 
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where  n  is  the  complex  index  of  refraction  defined  by(Budden,  1965) 


n  = 


2  1  2 
2  £  e2  sin  0  +  —  £3  !£1+e2)  (1+cos  9) 

(£l+£2>  Sin 
{sin4  8  £,£„  -  t 


2  2 
(e^+e^  sin  0  +  2  £3  cos  0 
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2  £3 


(£,+£n)  +COSQ  e3  (£!~E2)}  (2.106) 


*1  2 


2  2 
(£  +£  )  sin  8  +  2  £^  cos  9 


Equation  (2.106)  is  the  generalized  magnetoionic  formula  for  a  homo¬ 
geneous,  anisotropic  plasma  which  was  first  developed  by  Sen  and 

2 

Wyller  (1960) .  It  ls  apparent  that  r|  has  two  values:  one  correspond¬ 
ing  to  the  plus  sign  in  (2.106)  and  the  other  to  the  minus  sign.  The 
ionospheric  is  thus  a  birefringent  medium  which  supports  the  propa¬ 
gation  of  waves  of  two  different  types.  These  two  wave  types  corre¬ 
spond  to  two  different  polarizations  (e.g.  Budden,  1966)  and  are 
denoted  the  ordinary  or  0-mode,  HQ  (plus  sign) ,  and  the  extraordinary 


or  X-mode,  n  (minus  sign) . 

X 
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CHAPTER  III 
IONOSPHERIC  HEATING 
3 . 1  The  Energy  Balance  Equation 

The  objectives  of  this  chapter  are  to  determine  the  magnitude  and 
time  scale  of  electron  temperature  modifications  in  the  lower  iono¬ 
sphere  induced  by  a  strong  electromagnetic  wave  and  the  corresponding 
variation  in  the  heating  wave  fields  due  to  self -absorption. 

Electron  temperature  modifications  due  to  radiowave  heating  are 
governed  by  the  energy  balance  equation  (2.80) .  In  order  to  simplify 
the  analysis  of  this  equation,  it  shall  be  assumed  throughout  this 
chapter  that  the  time  scale  for  electron  density  modifications  is 
much  larger  than  that  for  the  electron  temperature  so  that  one  may 
write  Sn^/Bt  =  0.  The  validity  of  this  assumption  will  be  examined 
in  Chapter  IV.  Applying  this  assumption  to  (2.80),  the  energy  balance 
equation  becomes 

3  3  r  9 

-  n  K  -t—  T  =  E*J->3X-U  (3.1) 

2  e  B  3t  e  L  3t  es 

s 

where  E  •  J  is  the  rate  of  electron  heating  due  to  the  wave  field, 

and  £  3u  /3t  =  L  co11  is  the  rate  of  electron  energy  loss  due  to  col- 
s  es  e 

lisions  with  neutral  particles. 

From  equations  (2.16)  and  (2.88),  one  has 


E  •  J  =  E  •  (V  x  H  -  jo)£  E) 
e  o 


(3.2) 


=  -^  •  (E  x  H)  -  ju)(e  E2  +  y  H2) 

o  o 


By  computing  the  time  average  of  E  •  J  over  co/2tt,  one  finds  that  (e.g. 

e 


Brandstatter ,  1965) 


0 


(3.3) 


where 


S  =  R  {E  x  H*} 
2  e 


(3.4) 


is  the  average  power  density  of  the  heating  wave.  The  final  form  of 
the  production  term  (3.3)  to  be  used  in  the  solution  of  (3.1)  will  be 
developed  in  the  following  section. 

The  electron  energy  loss  term  )  9u  /9t  includes  a  number  of  col- 

s  es 

lisional  processes  including  elastic  electron-neutral  collisions  and 
electron  impact  excitation  of  the  rotational,  vibrational  and  electronic 
energy  levels  of  neutral  particles.  The  rate  of  electron  cooling  due 
to  elastic  electron-neutral  collisions  has  been  calculated  by  Banks 
(1966)  and  is  given  by 


J.  U  (T  )  |  _  _  “ 

3t  e  e  elastic  m 


n  /  V  (T  )\ 
e  \  es  e  / 


(T 
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T  )  . 
n 


(3.5) 


Elastic  collisions  are  not  very  effective  in  cooling  electrons  due  to 

the  unfavorable  mass  ratio,  m  /m  ,  and  these  processes  are  generally 

e  s 

overshadowed  by  inelastic  collisions.  At  present,  a  precise  analytic 

theory  of  inelastic  cooling  does  not  exist;  however,  numerical  formulae 

for  most  of  the  important  ionospheric  electron  cooling  processes  based 

on  laboratory  data  are  available,  and  have  been  summarized  in  works  by 

Stubbe  and  Varnum  (1972) ,  Banks  and  Kockarts  (1973)  and  Shunk  and  Nagy 

(1978) .  The  primary  electron  cooling  processes  in  the  ionosphere  for 

T  <  10  T  are  rotational  and  vibrational  excitation  of  N_  and  0.  and 
e  —  n  2  2 

the  excitation  of  fine  structure  levels  on  atomic  oxygen. 

Figures  3.1-2  illustrate  the  electron  temperature  variation  of 


T  =200 °K 
n 

1  -  ROTATIONAL  (Ng  ) 

2  -  ROTATIONAL  (02  ) 

3  -  ELASTIC  (N2) 

4  -  ELASTIC  (02  ) 

5  -  VIBRATIONAL  (N2  ) 
6-  VIBRATIONAL  (02  ) 


I/. 


8Te/Tn 

Figure  3.1  Normalized  electron  cooling  rates  for  various 
processes  in  air  as  a  function  of  the  frac¬ 
tional  change  in  electron  temperature. 
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the  electron  energy  transfer  rates  for  various  types  of  collisional 
processes  based  on  the  summary  of  formulae  given  by  Stubbe  and  Varnum 
(1972).  In  Figure  3.1,  the  mixing  ratio  for  the  major  neutral  gases 
has  been  assumed  to  be  constant  with  being  78%  of  the  total  neutral 
(air)  content  and  0^  being  21%.  The  cooling  rates  have  been  made  inde¬ 
pendent  of  the  electron  and  neutral  particle  densities  by  dividing 

9u  /3t  for  each  process  by  n  and  n  .  .  These  normalized  rates  are 
e  e  air 

plotted  as  a  function  of  the  normalized  change  in  electron  temperature 

above  ambient:  6T  /T  =  (T  -  T  )/T  .  For  6T  /T  £  1 ,  rotational 
e  n  e  n  n  e  n 

cooling  is  dominant;  however,  for  higher  electron  temperatures,  vibra¬ 
tional  cooling  becomes  important. 

In  the  D-region,  one  has  n  <<n^  and  n  ,  hence,  electron  cooling 

0  N2 

due  to  collisions  with  atomic  oxygen  is  relatively  unimportant,  but  in 

the  E-region,  n^  ~  n^  and  these  processes  must  be  considered.  Figure 
O  O2 

3.2  shows  a  plot  of  the  normalized  cooling  rates  for  electron  col¬ 
lisions  with  atomic  oxygen  as  a  function  of  <5t  /T  .  Electron  cooling 

e  n 

due  to  electron  impact  induced  fine  structure  transitions  in  0  is 
clearly  dominant  over  the  corresponding  elastic  cooling  process  over 
the  range  of  electron  temperatures  of  interest  in  this  study. 

3.2  The  Heating  Wave  Fields 

The  source  of  the  heating  wave  is  assumed  to  be  a  powerful  ground- 
based  HF  transmitter  which  is  excited  by  a  sinusoidally  varying  cur¬ 
rent  of  angular  frequency  u.  At  a  point  r  from  the  HF  transmitter 

(i.e.  r»l/K  )  but  below  the  base  of  the  ionosphere  (~  50  km)  the 
o 

heating  wave  fields  are  given  by  (Jordan  and  Balmain,  1968) 
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E  (r,  t)  =  -  j  K 


j  (U)t-K  r) 

z  £■ 

o  r 


I 


JEXT(r  ) 


-jK 


dr  (3.6) 


and  H  (r,  t)  =  \  r  x  E  (r,  t)  (3.7) 

o 

where  Z  =  /u  /e  -  characteristic  impedance  of  free  space, 
o  o  o 

and 

A 

r  -  a  unit  vector  in  the  direction  of  r. 

The  fields  incident  upon  the  ionosphere  therefore  have  the  form  of 
spherical  waves. 

In  deriving  the  magnetoionic  formula  (2.106) ,  it  is  assumed  that 
the  wave  fields  have  the  form  of  plane  waves.  If  one  were  instead 
to  assume  that  the  wave  fields  are  spherical,  then  the  application  of 
the  del  operator,  to  the  electric  field  in  (2.94)  gives  rise  to 
terms  of  order  1/r  as  well  as  those  of  order  K.  Clearly  if  one  has 
l/r<<K,  then  the  dispersion  equation  (2.96)  and  its  solution  (2.106) 
are  also  valid  for  spherical  waves.  At  the  base  of  the  ionosphere 
this  condition  reduces  to  co>>6  x  lO'Vs  which  is  easily  satisfied  at 
frequencies  in  the  HF  band.  Thus,  while  the  wave  fronts  of  the  radia¬ 
tion  field  from  the  transmitter  are  spherical,  by  the  time  they  reach 
the  base  of  the  ionosphere  the  radius  of  curvature  of  the  fronts  are 
so  large  that  the  wave  is  very  nearly  a  plane  wave.  The  wave  fields 
in  a  plasma  far  from  the  wave  source  are  therefore  of  the  forms 

j t cut  -  K  •  (r  -  r  )]  (3.8) 

e  o 

and 


E  (r,  t)  =  Eq  — 


H  (r,  t)  =  H  —  e 

O 


j  [u)t  -  K  •  (r  -  r  )] 


(3.9) 


where  E  is  the  value  of  E  at  some  reference  point  r  and  by  (2.16) 
o  o 
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Ho  up 


K  X  E 


(3.10) 


The  above  formulae  apply  to  a  homogeneous  plasma.  However, 
as  noted  in  Chapter  I,  the  ionosphere  is  a  nonhomogeneous  medium. 

For  a  nonhomogeneous  medium  the  spherical  wave  solutions  for  the 
field  vectors  given  above  are  not  valid  and  one  must  generally  resort 
to  a  full  wave  analysis  of  Maxwell's  equations  to  determine  the 
field  vectors  (e.g.  Budden,  1965)  .  At  wave  frequencies  above  1  MHz, 
the  properties  of  the  ionosphere  are  slowly  varying  over  distances  of 
the  order  of  a  wavelength.  One  may  therefore  consider  the  plasma  to 
be  locally  homogeneous  and  approximate  solutions  for  the  wave  fields 
can  be  obtained  by  the  WKB  method.  Since  the  largest  variation  in 
the  properties  of  the  ionosphere  is  in  the  vertical  direction,  these 
properties  are  assumed  to  be  a  function  of  height,  h,  only.  The  WKB 
solutions  for  the  portion  of  the  spherical  wave  front  that  is  verti¬ 
cally  incident  on  a  horizontally  homogeneous  ionosphere  are  given 
approximately  by  (Budden,  1966) 


E  (h,  t)  =  E 


-h, 


and 


n  (h)  exp  /  jut  -  j  K o  j 


H  (h,  t)  =  Hq  I  ^  l")*5  (h)  exp  |  jut  -  j  Ko 


n 

! 


n(h')  dh' 


(h  )  dh 


(3.11) 


(3.12) 


where 


H  =  z  x  E  /Z 
o  o  o 


(3.13) 


Substituting  (3.11-13)  into  (3.4),  one  finds  that  the  average 
power  density  is  given  approximately  by 
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S  (h)  =  S 


o  \  h 


n 

r 


exp  (  -2Kq  |  x(h  )  dh 


(3.14) 


where  X(h)  is  the  imaginary  part  of  the  index  of  refraction 


n (h)  =  U(h)  -  j  x(h) 


(3.15) 


and  u(h)  is  the  real  part.  The  quantity  Sq  is  the  average  power  den¬ 
sity  which  is  vertically  incident  on  the  base  of  the  ionosphere,  h  , 
and  is  given  by 


S 

o 


(3.16) 


where  PT  -  total  transmitter  power,  and 

G  -  zenith  gain  of  the  transmitting  antenna  array. 
Taking  the  divergence  of  equation  (3.14)  ,  one  finds  that 


^  •  S(h)  =  "rS(h)  -  2  K  X(h)  S(h) 
h  o 


(3.17) 


The  first  term  on  the  righthand  side  of  (3.17)  accounts  for  the  decrease 
in  the  power  density  of  the  HF  wave  as  its  front  spreads  with  increas¬ 
ing  height,  while  the  second  term  represents  the  power  dissipation  per 
unit  volume  as  a  result  of  absorption  by  electrons.  The  rate  of  elec¬ 
tron  heating  due  to  the  wave  field  is  therefore  given  by 

E  •  J  =  2  K  x (h)  S  (h)  (3.18) 

e  o 

Substituting  this  result  into  the  electron  energy  balance  equation, 
one  obtains 


4  n  (h)  K  37  T  (h)  =  2  K  X  T  )  S  (h,  T  ) 
2e  Bdte  o  e  e 


-  I  37  0  (h,  T  ) 

L  dt  e,s  e 

s 


L 


(3.19) 
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where  the  electron  temperature  dependence  of  x  and  s  arises  from  the 

collision  frequency,  v  (K  T  ) ,  which  is  a  parameter  in  the  formula 

en  B  e 

for  the  index  of  refraction. 

As  a  heating  wave  propagates  through  the  ionosphere,  energy  is 
absorbed  from  it  by  the  plasma  electrons,  causing  an  increase  in  the 
electron  temperature  as  dictated  by  (3.19).  The  increase  in  electron 
temperature  changes  the  plasma  absorption  index  from  its  ambient  value, 
X(h,  T  )  to  a  heated  value,  T  )  •  along  the  wave  path  which  in 

turn  causes  additional  absorption  on  the  wave.  From  (3.14),  one  finds 
that  the  self-absorption  experienced  by  a  heating  wave  as  it  propagates 
from  hQ  to  h  is  given  by 

h 

Ls  =  8.686  Kq  J  [X(h",  Tg)  -  X(h^  TM]  dh'  (3.20) 

h 

O 

where  Ls  is  in  decibels  above  the  ambient  absorption  level. 


3 . 3  Electron  Heating  in  a  Weak  Field:  Analytic  Theory 

In  the  following  analysis,  the  heating  wave  field  is  assumed  to 

be  weak  so  that  the  electron  temperature  is  only  slightly  perturbed 

from  its  ambient  value:  6  T  <<  T  .  In  this  temperature  regime,  the 

e  n 

dominant  electron  energy  loss  process  is  rotational  cooling  with  N2 
and  C>2  (see  Figure  3.1),  and  the  total  cooling  rate  can  be  written  as 
(Mentzoni  and  Row,  1963;  Shunk  and  Nagy,  1978) 


7  ^-U  (T  )  =  a  n  T  6t 

*•  at  e,s  e  e  e  e 

s 


(3.21) 


where 


a  =  1.17  x  10-10  (.78  n,  of  B  N2  +  .21  n,  qI  B  °2) 

“2  ^2  °  °2  °2  ° 
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Since  T  =  T  ,  one  can  neglect  self -absorption  effects  and  the  electron 
e  n 


energy  balance  equation  can  be  written  as 

.  St 

3  T  =  P  (T  )  -  ® 


where 


and 


3t  e  e  '  n  t  (T  ) 
ri  n 


4  K  X<T  )  S (T  ) 

,  on  n 

P  (T  )  =  - r— — - 

e  n  3  n  K_ 


e  B 


3  K  T  ' 
,  Bn 

TH  Tn  = 


(3.22) 


(3.23) 


(3.24) 


For  an  HF  wave  of  duration  At  initiated  at  time  t  =  0,  one  finds 

H 

from  equation  (3.22)  that  for  0  <  t  <  At 

H 

T  (t)  =  T  +  P  (T  )  T  (T  )  [1  -  e"t/TH(Tn)] 


n  e  n  H  n 


where 


St  (At  )  =  P  (T  )  T  (T  )  [l  -  e'Atli/TH  {Tn)  ] 
e  H  e  n  H  n 


(3.25) 


(3.26) 


3.4  Electron  Heating  in  a  Strong  Field:  Computational  Method 

This  section  describes  a  numerical  procedure  for  solving  the 
energy  balance  equation  for  a  strong  wave  field  in  a  nonhomogeneous 
plasma  which  takes  into  account  the  self-absorption  of  the  heating 
wave. 

Let  h  denote  the  base  of  the  ionosphere  and  let  the  plasma  above 
o 

this  height  be  subdivided  into  a  large  number  of  horizontal  layers  of 
thickness  Ah.  =  h.  ,  -  h.  where  h..,  and  h.  are  the  height  of  the 
upper  and  lower  boundaries  of  the  i-th  layer.  Each  Ah^  is  chosen  such 
that  the  plasma  is  essentially  homogeneous  over  the  layer.  In  partic¬ 
ular,  it  is  required  that  Ah.  <<  1/(2  K  Ay.)  where  Ay.  =  Y(h.  ,)-Y(h.) 

1  0^1  !  i+l  A  i 
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so  that  the  incremental  absorption  over  each  h^  is  small  compared  to 

the  integral  absorption  from  hQ  to  h^. 

Consider  an  HF  pulse  of  duration  At  incident  upon  h  at  t 

H  o  o 

=  hQ/c  and  propagating  upwards  with  velocity  c.  If  one  chooses 

t.  «  r  (h.  ,  T  )  where  At.  is  the  size  of  the  time  increments  within 
1  H  1  n  l 

the  i-th  layer,  the  electron  temperature  within  the  i-th  layer  can  be 
approximated  as 


T  (h. ,  t.)  =  T  (h. , 
e  i  3  e  i 

where  t .  =  t .  ,  +  At . , 
3  3-1  i 


t.  >  +  At.  fie  (h.,  t.  ,,  S(h.,  t.))  (3.27) 
3-1  l  it  1  3-1  1  3 

and  3Tg/3t  is  given  by  equation  (3.19)  with  X 


and  7  3u  /3t  evaluated  at  T  (h.,  t.  ,).  T  is  computed  for  each  t. 
g  es  e  l  3-1  e  3 

starting  at  hQ  and  moving  upwards  in  steps  of  Ah^.  The  power  density 

is  recomputed  after  each  altitude  increment  for  fixed  t_.  using  equation 

(3.14)  and  is  employed  in  evaluating  Tg  for  the  layer.  The  boundary 

conditions  on  the  solutions  for  T  and  S  are 

e 


S (h  ,  t) 
o 


S  ,  t  <  t  <  t  +  Atu 
o  o  —  —  o  H 

0,  otherwise 


and 


(3.28) 


T  (h.,  t)  =  T  (h.),  t  <  h./c  (3.29) 

e  i  n  l  l 


3.5  Electron  Heating  in  a  Strong  Field:  Results 

The  above  procedure  has  been  used  to  compute  the  transient  and 

steady  state  response  of  T  and  S  for  a  strong  wave  field.  The  absorp- 

e 

tion  index  was  calculated  from  the  Sen-Wyller  formula  (2.106) .  The 

-► 

gyro  frequency  was  taken  as  f  =  1.135  MHz  and  the  angle  between  K 

G 

and  Y,  0  =  40.65°,  corresponding  to  geomagnetic  conditions  at  Arecibo 


for  a  vertically  incident  wave.  The  collision  frequency  was  calculated 

from  the  energy  dependent  cross  sections  given  by  Banks  and  Kockarts 

(1973)  and  was  evaluated  at  the  most  probably  electron  energy,  K^T^. 

The  neutral  particle  densities  and  the  neutral  temperature  were  taken 

from  CIRA  (1972) .  The  plasma  frequency  model,  U)  (h) ,  corresponds  to 

e 

the  daytime  electron  density  distribution  at  a  solar  zenith  angle,  Xg 
=  30°  and  a  moderate  level  of  solar  activity. 

Figure  3.3  shows  the  electron  temperature  variation  at  70  km  due 
to  a  250ys,  5  MHz  X-mode  heating  pulse  for  several  values  of  the  effect¬ 
ive  radiated  power,  ?TG.  One  notes  that  the  electron  temperature 
reaches  a  steady  state  in  less  than  250ys  for  all  heating  powers  and 
that  the  plasma  cools  somewhat  faster  than  it  heats.  The  power  density 
at  70  km  (not  shown)  is  proportional  to  the  heating  power  and  is  con¬ 
stant  with  time.  This  is  not  the  case  in  the  upper  D-region  (90  km) 
as  illustrated  in  Figure  3 . 4.  Here ,  the  power  density  drops  abruptly 
from  its  initial  value,  corresponding  to  the  leading  edge  of  the  heating 
pulse,  to  a  steady  state  value  on  a  time  scale  of  about  500ys.  This 
effect  is  due  to  self -absorption  and  it  parallels  the  time  scale  for 
electron  temperature  changes  in  the  70  to  80  km  region,  which  is  the 
region  of  maximum  heating  as  illustrated  in  Figure  3.5.  The  electron 
temperature  at  90  km,  on  the  other  hand,  responds  much  more  slowly 
than  the  power  density,  requiring  about  3.5  ms  to  reach  a  steady  state 
level . 

Self-absorption  severely  limits  the  available  heating  power  in  the 
upper  D-region  leading  to  a  steep  temperature  gradient  between  80  and 
90  km  as  shown  in  Figure  3.5.  Moreover,  the  limitations  imposed  by 


Power  density  and  electron  temperature  variation  at  90  km  due  to 
heating  pulse. 
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self-absorption  on  upper  ionospheric  heating  cannot  be  readily  over¬ 
come  simply  by  increasing  the  heating  power.  In  Figure  3.6,  the 
initial  (leading  edge)  and  steady  state  power  densities  are  plotted 
as  a  function  of  heating  power  at  5  MHz  x-mode.  The  initial  power 
density  increases  in  proportion  to  the  heating  power >  however,  the 
self-absorption  also  increases  with  heating  power,  with  the  effect 
that  the  net  gain  in  steady  state  power  density  is  very  small.  For 
example,  at  100  km  the  self-absorption  is  about  3dB  for  100  MW  of 
heating.  If  one  doubles  the  heating  power  to  200  MW,  the  self¬ 
absorption  rises  to  more  than  5dB  resulting  in  a  net  gain  in  power 
density  of  less  than  30%. 

Figure  3 . 7  shows  the  initial  and  steady  state  power  density  as 
a  function  of  heating  frequency.  The  power  density  available  for 
heating  in  the  upper  ionosphere  decreases  sharply  at  the  lower  end  of 
the  heating  frequency  range  due  to  rapid  absorption  in  the  middle  D- 
region.  As  a  consequence,  the  electron  temperature  in  the  middle  D- 
region  is  greatly  enhanced  (Figure  3.8).  As  one  moves  to  higher 
heating  frequencies,  the  D-region  absorption  significantly  decreases 
and  consequently  so  does  the  D-region  electron  temperature.  The  self¬ 
absorption  component  of  the  D-region  absorption  also  decreases  with 
increasing  heating  frequency  as  shown  by  the  convergence  of  the  initial 
and  steady  state  power  density  curves  in  Figure  3.7.  At  the  higher 
frequencies,  the  power  density  available  for  upper  ionospheric  heating 
becomes  nearly  independent  of  frequency  and  is  controlled  by  the  heating 
power.  While  the  available  power  density  is  much  larger  at  the  higher 
frequencies,  the  heating  above  90  km  is  not  significantly  increased,  as 
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shown  in  Figure  3.8,  because  the  critical  power  density,  S  ,  required 

c 

to  produce  large  temperature  enhancements  increases  with  heating  fre- 

2 

quency,  varying  approximately  as  f  (Meltz  et  al.,  1974). 

The  results  presented  above  can  be  influenced  by  a  number  of 

factors.  The  effect  of  two  of  these  factors,  variations  in  electron 

density  and  uncertainties  in  the  total  electron  energy  loss  rate,  are 

considered  in  Figure  3.9.  The  ambient  (T  =  T  )  temperature  model  and 

e  n 

that  for  100  MW  ERP ,  5  MHz  X-mode  are  illustrated  as  curves  (1)  and 

(2)  ,  respectively.  Curves  (3)  and  (4)  show  the  effect  of  a  factor  of 
two  changes  in  the  electron  density  at  all  heights  on  the  electron  tem¬ 
perature  during  heating.  Increasing  the  total  electron  concentration 

(3)  results  in  only  a  moderate  decrease  in  T  ,  primarily  near  the 
temperature  maximum  (73  km) .  On  the  other  hand,  decreasing  n^  (4) 
causes  a  considerable  enhancement  in  electron  temperature  in  the  upper 
D-region,  which  suggests  that  the  upper  ionosphere  is  heated  more  ef¬ 
fectively  as  the  solar  zenith  angle  increases.  Curves  (5)  and  (6) 

illustrate  the  effect  of  a  factor  of  2  increase  and  decrease  in  the 

v*  9 

total  electron  energy  loss  rate,  )  q—  U  ,  on  the  electron  temperature 

§  dt  e , s 

distribution.  The  electron  temperature  maximum  varies  roughly  as  the 

inverse  of  the  total  loss  rate.  This  result  could  be  anticipated  from 

the  analysis  given  in  Section  3.3.  From  (3.25)  one  expects  the  steady 

state  electron  temperature,  T  ,  for  the  strong  field  case  to  be  given 

eo 

approximately  by 

T  =  P  (T  )  T  (T  ) 
eo  e  eo  H  eo 


(3.30) 


IGHT(Km) 


ELECTRON  TEMPERATURE  (K) 

Figure  3.9  Effect  of  electron  density  and  total  loss 

rate  on  the  electron  temperature  distribution. 
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where  the  analogous  form  of  (3.24)  is 


r  (T  )  =  ~  n  K  St  /  T  3u  /9t 
H  eo  2  e  B  eo  5  as 


(3.31) 


Thus,  T  is  inversely  proportional  to  the  total  electron  energy 


loss  rate. 


CHAPTER  IV 


ION  CHEMISTRY  MODIFICATIONS 
4.1  Continuity  Equations 

In  the  previous  chapter,  it  was  shown  that  a  substantial  increase 
in  electron  energy  can  occur  when  the  ionosphere  is  irradiated  by  high 
power  electromagnetic  waves.  Since  the  chemical  processes  which  deter¬ 
mine  the  composition  of  the  ionospheric  plasma,  and  in  particular,  the 
electron  density,  are  collisional  in  nature,  they  depend  directly  upon 
the  energies  of  the  colliding  particles.  An  increase  in  electron  energy 
due  to  heating  by  electromagnetic  waves  can  therefore  alter  the  reaction 
rates  for  chemical  processes  in  which  electrons  participate  and  thus 
modify  the  composition  of  the  plasma. 

In  Chapter  III  it  was  assumed  that  the  time  scale  for  HF  heating 

induced  electron  density  variations  is  much  larger  than  that  for  the 

electron  temperature  so  that  one  could  write  Bn  /Bt  =  0  in  the  electron 

e 

energy  balance  equation.  The  objectives  of  this  chapter  are  to  ascertain 
the  validity «of  this  assumption  and  to  explore  the  types  of  chemistry 
modifications  that  can  result  from  HF  heating. 

The  electron  content  of  the  ionospheric  plasma  is  governed  by  the 
continuity  equation  (2.82)  which  can  be  written  out  explicitly  in  terms 
of  the  various  electron  production  and  loss  processes  as 


+ 


l 


Y.  n.  -  6  (T  )  n 
i-  1-  e  e  e 


l 

j+ 


atJ_  (T  ) 
3+  e 


n . 

3  + 


n 

e 


(4.1) 


Q 


e 


where 


-  electron  production  rate  per  unit  volume  due  to 
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ionization  of  neutral  atmospheric  constituents, 

Yi  -  rate  of  electron  detachment  from  the  i-th  negative 

ion  type , 

8  (T  )  -  rate  of  electron  attachment  to  neutral  constituents, 

e  e 

and 

a.  (T  )  -  rate  coefficient  for  electron  recombination  with 

3+  e 

the  j-th  positive  ion  type. 

The  electron  attachment  rate  and  the  electron-ion  recombination  coef¬ 
ficients  depend  directly  on  the  electron  temperature  as  emphasized  in 
(4.1) .  Strictly  speaking,  the  electron  production  rate  may  also  depend 
on  T^  due  to  electron  impact  ionization  of  neutral  particles  (Gurevich 
et  al.,  1977b),  or  may  depend  directly  on  the  electric  field  strength 
of  the  heating  wave,  E.  The  latter  dependence  results  from  the  break¬ 
down  of  the  neutral  gas  in  the  presence  of  very  strong  electric  fields 
(Lombardini,  1965) .  The  effects  of  electron  production  enhancements 
due  to  the  heating  wave  field  shall  be  commented  on  later  in  this  chap¬ 
ter.  In  the  present  discussion,  it  is  assumed  that  Q  is  independent  of 

e 

T^  and  E  and  is  given  by  the  sum  of  the  production  rates  due  to  ioniza¬ 
tion  of  the  neutral  gas  by  solar  photons  and  cosmic  rays. 

In  the  upper  D  and  lower  E  regions,  the  predominant  positive  ions 
are  NO  and  O^*,  there  are  very  few  negative  ions,  and  electron  attach¬ 
ment  and  detachment  processes  are  relatively  unimportant.  From  the 
plasma  charge  neutrality  condition,  p  =  0,  one  has 


n 

e 


NO"' 


+  n 


°2+ 


(4.2) 
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In  this  region  the  electron  continuity  equation  can  be  written  in  the 
simple  form 


_3 

at 


n 

e 


a 


eff 


(T  )  n  2 
e  e 


(4.3) 


where  a  ,^(T  )  is  the  effective  electron-ion  recombination  coefficient 
eff  e 

defined  by 


a  (T  ) 
eff  e 


=  I 

j  + 


a.  AT  ) 
3+  e 


V? 

J  1+ 


ni+ 


(4.4) 


Holway  and  Meltz  (1973)  have  employed  an  analytic  model  of  this  form  to 
calculate  the  electron  density  variation  due  to  various  levels  of 
heating.  They  have  found  that  the  electron  density  increases  from  its 
unheated  value  on  a  time  scale  of  minutes  to  a  level  which  is  inversely 
proportional  to  the  square  root  of  the  recombination  coefficient. 

While  simple  models  of  plasma  density  modification  based  strictly 

upon  heating  induced  changed  in  the  recombination  rate  may  be  adequate 

in  the  upper  D  and  lower  E  regions,  the  complex  nature  of  the  ion  chem- 

istry  governing  the  plasma  composition  in  the  lower  D-region  requires 

a  more  detailed  analysis.  Below  85  km  the  identity  of  the  predominant 

positive  ions  switches  from  N0+  and  02+  to  water  cluster  ions  of  the 

forms  H+* (HO)  and  N0+# (HO)  which  are  formed  through  a  complicated 
2  n  2  m 

sequence  of  reactions  (e.g.  Rowe  et  al.,  1974).  Electron  recombination 
with  these  heavier  ions  will  also  be  influenced  by  HF  heating.  More¬ 
over,  below  75  km  negative  ions  become  a  significant  portion  of  the  total 
ion  density  due  to  the  attachment  of  electrons  to  neutral  constituents, 
primarily  through  the  three  body  reaction  e  +  0^  +  M  -*■  0^  +  M.  The 


rate  for  this  reaction,  $,  should  increase  with  electron  temperature 
up  to  about  700°K  for  ambient  neutral  temperatures  of  200-300°K  then 
slowly  decrease  (Chanin  et  al.,  1962).  Thus,  HF  heating  of  the  lower 


D-region  could  decrease  the  electron  density  due  to  increased  three 
body  attachment.  However,  this  process  must  compete  with  any  electron 
increase  due  to  decreased  recombination. 

Equation  (4.1)  is  inadequate  for  describing  the  chemically  rich 
lower  ionosphere  since  some  of  the  electron  production  and  loss  terms 
depend  on  the  ion  densities  which  are  unknown.  One  must  therefore 
develop  a  set  of  continuity  equations  for  the  ionic  components  of  the 
plasma.  The  ion  continuity  equations  can  be  derived  directly  from 
Boltzmann's  equation  (2.20)  by  integrating  over  dv^  (e.g.  Holt  and 
Haskell,  1965) .  Neglecting  transport  effects,  one  finds  that  the  i-th 
type  of  positive  ion  is  described  by  a  continuity  equation  of  the  form 


ni+  =  +  I  (R.  ,  n  -  R .  n.  ) 

dc  i+  i+,3+  3+  3+,i+  i+ 


(4.5) 


-(oc.(T)n  +  I  a  .  n.  )  n. 


i+  e  e  i+,3~  3-  1+ 


while  the  i-th  type  of  negative  ion  is  governed  by  a  continuity  equation 
of  the  form 


8  n . 


-57-  “i-  =  6.  (T  )n  +  T  (R.  n.  -  R.  n.  ) 

dt  1-  e  e  i-,3-  3~  3~,i-  1- 
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where 

Q^+  -  production  rate  per  unit  volume  of  the  i-th  positive 

ion  type  resulting  from  the  ionization  of  neutral 
atmospheric  constituents  by  solar  radiation  and 
cosmic  rays, 

R.+  -  rate  at  which  the  i-th  positive/negative  ion  type  is 

created  from  the  j-th  positive/negative  ion  type 
through  charge  transfer  and  atom-ion  interchange 
with  neutral  constituents, 

a.  -  ion-ion  recombination  rate  coefficient  for  the  i-th 
r+fD- 

positive  ion  type  and  j-th  negative  ion  type,  and 

8^  -  rate  at  which  the  i-th  negative  ion  type  is  formed 

by  the  attachment  of  electrons  to  neutral  consti¬ 
tuents  • 

The  system  of  equations  (4.1) ,  (4.5-6)  are  not  independent,  but  are 
related  by  the  plasma  neutrality  condition 


n  =  £  n .  -  £  n. 

e  ",  1+  r  i- 


(4.7) 


l- 


Equations  (4.5-7)  specify  the  composition  of  the  plasma  provided  one 
knows  the  electron  temperature  distribution.  As  shown  in  Chapters  IX  and 
III,  the  electron  temperature  distribution  depends  on  the  heating  wave 
field  strength,  which  in  turn  depends  on  the  electron  temperature  and 
electron  density  through  the  index  of  refraction  of  the  plasma.  Therefore, 
a  complete  description  of  heating  wave  induced  plasma  modification 
requires  a  simultaneous  solution  of  the  system  of  equations  (4.5-7) ,  the 
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electron  energy  balance  equation,  and  the  heating  wave  power  density 
equation  (see  Chapter  III)  .  A  simultaneous  solution  of  this  set  of 
equations  is  complicated  by  the  fact  that  the  lower  ionosphere  con¬ 
tains  a  large  number  of  different  types  of  ions  and  many  of  the  reac¬ 
tion  rates  leading  to  the  formation  and  destruction  of  these  ions  are 
not  known.  Accordingly,  the  following  section  describes  a  simplified 
ion  chemistry  scheme  which  avoids  some  of  the  complexity  and  uncertainty 
inherent  in  detailed  models  of  the  lower  ionosphere. 

4.2  The  Mitra-Rowe  Ion  Chemistry  Model 

The  ion  chemistry  model  adopted  in  this  work  is  the  simple  six  ion 
scheme  proposed  by  Mitra  and  Rowe  (1972) .  This  model  has  been  shown  to 
satisfactorily  reproduce  the  major  features  of  the  lower  ionosphere 
during  quiet  and  disturbed  conditions  while  eliminating  many  of  the 
complex  and  often  hypothetical  steps  in  the  complete  ion  chemistry 
reaction  sequence  (Mitra,  1975a,  b;  Rowe  et  al.,  1974).  The  reaction 
chain  for  the  six  ion  model  is  outlined  in  Figure  4.1  and  specific 
reactions  included  in  the  model  are  listed  in  Table  4.1. 

The  primary  ions  formed  by  the  ionization  of  neutral  atmospheric 

constituents  are  N2+,  02+  and  N0+.  The  ions  N2+  and  C>2+  are  formed  by 

two  different  mechanisms :  photoionization  of  the  neutral  gas  by  solar 

X-rays  and  extreme  ultraviolet  radiation,  and  collisional  ionization 

of  N2  and  02  by  high  energy  cosmic  particles  (i.e.  cosmic  rays) .  Nitric 

oxide  is  a  minor  neutral  constituent  (i.e.  n  _  «  n_  or  n ,  )  but  is 

NO  O2 

very  effectively  ionized  by  the  solar  Lyman  a  line  at  1215.7  A,  forming 


NO  . 
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Table  4.1  Reactions  included  in  the  six  ion  scheme 


Ki 


1.  N^+  +  "*■  02+  +  NjJ  =  5.0  x  10  ^  cm3/s 

2.  N2+  +  0  -*•  NO+  +  N;  K  =  1.4  x  10_1°  cm3/s 


+ 

+ 

-10 

1  3. 

°2 

+  NO 

-►  NO 

+  °2  ' 

K3  =  4.4  x 

10 

+ 

+ 

,  -16 

3 

1  4. 

°2 

+  n2 

-►  NO 

+  NO; 

K.  =  10 

4 

cm  /s 

5.  0n+  +  0  +  0„  -»■  O*  +  0_;  K  =  2.8  x  10“3° (300°K/T  )cm3/s 

2  2  2  4  2  5  n 

6.  04+  +  O  -*■  02+  +  03;  Kg  =  3  x  10_1°  cm3/s 

7.  04+  +  H20  -*■  02  +  02+  •  H2Oi  K7  =  1.5  x  10  9  cm3/s 

8.  02  +  o  -*■  C>3  +  e;  Kg  =  3  x  10  10  cm3/s 

9-  02-  +  02  (■'■Ag)  -*■  2  02  +  e;  Kg  *  2  x  lo"10  cm3/s 

10.  hv  +  o  “  0,  +  e;  Y  =0.3  s_1 

2  2  p 


( 
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The  rate  of  production  of  the  i-th  ion  type  per  unit  volume  at 
height  h  due  to  the  photoionization  of  the  i-th  type  of  neutral  con¬ 
stituent  by  all  solar  radiations  is  given  by  (e.g.  Rishbeth  and 
Carriott,  1969;  Banks  and  Kockarts,  1973) 


(h)  =  a.  (X)  n.(h)  $(h,  XI  dX" 


where 

(X)  -  the  ionization  cross  section  at  wavelength  X* 

for  the  i-th  neutral  constituent,  and 
$(h,  X")  -  the  solar  radiation  flux  of  wavelength  X  a 
height  h. 

For  ionizing  collisions  of  the  i-th  neutral  particle  type  with  high 
energy  particles  (i.e.  >  100  KeV)  the  corresponding  production  rate  is 
given  by  (Velinov,  1968) 
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f  (u)  -  the  energy  spectrum  of  the  j-th  type  of  energetic 
particle  at  the  top  of  the  earth's  atmosphere. 

The  production  rates  for  N^"1",  ®2+  an<^  N®+  use<^  ^is  wor^  have  been 
calculated  from  the  above  formulae  using  a  computer  code  developed 
by  Rowe  (1972)  .  Typical  daytime  values  of  Q^+/  Qq+  and  0^+  are  given 
in  Figure  4.2. 

Since  N^+  is  rapidly  converted  to  0^+  or  N0+  by  reactions  (1)  and 
(2)  of  Table  4.1,  only  0 *  and  NO+  are  explicitly  included  in  the 
chemistry  model.  The  effective  production  rates  for  these  ions  are 


denoted  Qq+  and  Q^0+»  respectively,  and  are  defined  by 


and 


where 


202+  ”  202+  +  1  +  A  ^N2+ 
2NO+  “  %0+  +  1  +  A  2N2+ 

a  =  k2[o!/k1[o2] 


(4.10) 

(4.11) 


The  brackets  denote  the  number  density  of  the  indicated  constituent 
(e.g.  [o2]  =  nQ  ) . 

In  addition  to  the  primary  ions  and  N0+,  the  positive  ion 

scheme  includes  the  intermediate  ion  04+  and  a  hydrated  ion  group  Y+, 

which  includes  H+(H_0)  and  N0+(H„0)  .  The  primary  negative  ion  formed 
2  n  2  m 

by  electron  attachment  is  C>2-.  Al1  other  negative  ions  are  assumed  to 
be  formed  by  reactions  with  0^“  and  are  lumped  into  the  group  X  . 

The  rates  for  conversion  of  one  positive  ion  type  into  another  as 


indicated  in  Figure  4.1  are 


ION  PRODUCTION  RATE  (cm3  s'1) 

production  rates  as  a  function  of  height. 
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*3  t-°]  +  K4  [N2] 

(4.12) 

R2  =  K5  ^°2^2 

(4.13) 

R3  ‘  K6  f°> 

(4.14) 

'A  '  S  lH20!- 

(4.15) 

The  bulk  rate  for  conversion  of  N0+  into  H+ (h  O)  ,  B,  has  been  estimated 

2  n 

from  positive  ion  composition  measurements  and  is  given  approximately 
by 

B  =  10-31  [air]2  in  s-1  (4.16) 

(Rowe  et  al.,  1974).  Other  parameters  included  in  the  model  are  the 
electron  detachment  rate  from  , 

=  Yp  +  Ka  t°]  +  Kg  [°2(1Ag)],  (4.17) 

the  bulk  rate  for  conversion  of  0^~  into  X", 

X  =  10  [02 ]  [air]  +  3  x  10  ^  [o^ ]  in  s  3  (4.18) 

(Mitra,  1975b) .  The  ion-ion  recombination  coefficients  are  assumed  to 

be  the  same  for  all  ions  =  10  ^  cm3/s. 

One  addition  parameter  employed  in  the  model  is  the  bulk  detachment 

rate  from  X  ,  y^.  Rowe  et  al.  (1974)  have  shown  that  y^  plays  a  key 

role  in  determining  the  size  of  the  electron  density  and  negative  ion 

to  electron  ratio,  X,  below  70  km.  Estimates  of  y^  at  60  km  based 

upon  electron  and  total  positive  ion  density  measurements  range  from 

0.1s  3  to  4s  The  effect  of  y_  on  n  and  X  is  illustrated  in 

2  e 
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Figure  4.3.  The  solid  curves  represent  the  ambient  profiles  which 
will  be  adopted  in  subsequent  calculations.  The  corresponding  values 
of  y2  are  given  in  Table  4.2.  The  dashed  curves  in  Figure  4.3  are 
for  fixed  y^  as  indicated. 

Electron  attachment  in  the  lower  ionosphere  is  believed  to  take 

place  primarily  through  the  three  body  reaction  e  +  0^  +  -*■  0^  +  °2‘ 

A  large  number  of  experimental  studies  of  this  reaction  have  been 

performed,  and  it  is  perhaps  the  best  documented  of  all  negative  ion 

reactions  (Caledonia,  1975) .  For  unheated  conditions,  T  *  T  ,  the 

e  n 

attachment  coefficient  for  this  reaction  is  given  by  (Truby,  1972) 

K  (T  )  =  4  x  10-30  e“196/Tn  cm6/s;  100  <  T  <  600°K  (4.19) 

an  n 

Measurements  of  the  electron  energy  dependence  of  the  three  body  process 

for  fixed  values  of  T  (Chanin  et  al.,  1962;  Pack  and  Phelps,  1965)  have 

n 

shown  that  K  increases  rapidly  for  /u\  >3K  T  /2,  reaching  a  peak  at 
d  \  •  B  n 

an  average  electron  energy  of  (u^  =  0.09  eV  and  then  decreasing  slowly 
at  higher  energies.  The  0.09  eV  peak  corresponds  to  the  fourth  vibra¬ 
tional  level  of  02_*  resolution  electron  beam  measurements  by 

Spence  and  Schulz  (1972)  have  shown  that  has  secondary  maxima  at 
energies  corresponding  to  the  higher  vibrational  levels  of  C>2_. 

The  above  mentioned  experimental  data  has  been  used  to  construct 

a  numerical  formula  for  the  simultaneous  T  and  T  dependence  of  K 

e  n  a 

(see  Appendix  A) .  The  electron  temperature  variation  of  K  computed 
from  this  formula  is  compared  with  the  experimental  data  in  Figure  4.4, 
and  is  in  good  agreement.  The  electron  attachment  rate  employed  in  the 


Figure  4.3  Effect  of  y  on  the  electron  density  and  the  negative  ion  to  electron  ratio 
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Mitra-Rowe  chemistry  scheme  is  given  by 


8<T  ,  T  )  =  K  (T  ,  T  )  [0,]‘ 
e  n  a  e  n  i 


(4.20) 


Reviews  of  the  available  experimental  data  on  atmospheric  electron- 
ion  recombination  processes  have  been  given  by  Biondi  (1969,  1975) . 

In  the  lower  ionosphere  electron-ion  pairs  are  removed  by  the  highly 
efficient  dissociative  recombination  process.  Estimates  of  the  dis¬ 
sociative  recombination  coefficients  for  0^+  and  NO+  are 


where 


and 


where 


0t  +  (T  )  =  +  (T  )  (T  /T  ) 

0^  e  n  e  n 


-0.6 


a  .  (T  )  =  2.1  x  10~7  (T.  /300°K)"0*7  cm3s'1 
°2  n  n 


a  +  (T  )  =  Ot  .  (T  )  (T  /T  ) 
trcr  e  N0+  n  e  n 


-0.5 


(4.21) 


(4.22) 


(4.23) 


(T  )  =  4.1  x  10-7  (T  /300°K)  cm3s-1;  200  <.  T  <.  300°K  (4.24) 

no  n  n  n 


More  recent  measurements  of  dissociative  recombination  cross  sections 

for  0  +  and  NO+  by  Mul  and  McGowan  (1979)  suggest  that  both  a  +  and 

2  °2 
m  -0.5 

aNQ+  vary  as  Tg 

The  dissociative  recombination  coefficient  for  04+  is  approximately 

a  .  =  2.3  x  10  6  cm3s  1  at  T  =  T  =  205°K  and  is  assumed  to  vary  as 
°4  e  n 

T  • 0 . 
e 

The  effective  recombination  coefficient  for  the  water  cluster  ion 
group ,  Y+,  is  given  approximately  by 


88 


«v+  CTJ  =  I  (TJ 

n 


YT  "e'  *  n  -n'  [  Y+  ] 


f  H  •  (H  O)  I  T  -b 
1 _ 2  n  J  e 

T 


(4.25) 


where  a  (T  )  denotes  the  recombination  coefficient  for  the  hydrate 
n  n 

H+  •  (H_0)  as  given  by  Leu  et  al.  (1973),  with  a  neutral  temperature 
2  n 

-0 . 2 

dependence  of  '  .  The  electron  temperature  dependence  of  the 
effective  recombination  coefficient  is  uncertain,  but  typical  estimates 
are  0  <  b  <  1.  Recent  measurements  by  Huang  et  al.  (1978)  favor  a 
small  value  of  b  =  0.08.  The  relative  concentration  of  the  individual 
hydrated  ions  [H+  *  (H^O) q] /[Y+]  may  be  calculated  from  the  equilibrium 
rate  constants  for  the  reactions  H_0  +  H+  •  (H  o)  ,  H+  •  (H„0)  as 
given  by  Kebarle  et  al.  (1967) . 

The  reaction  rates  used  in  the  Mitra-Rowe  chemistry  scheme  vary 
as  a  function  of  altitude  due  to  their  dependence  on  the  neutral  particle 
densities.  The  neutral  atmospheric  model  adopted  in  this  work  is 
illustrated  in  Figure  4.5.  The  densities  of  the  major  neutral  consti¬ 
tuents  N2  and  C>2  are  taken  from  CIRA  (1972) ,  while  the  concentrations 
of  the  minor  constituents  0,  H20,  02(1Ag) ,  0^  and  NO  are  based  upon 
the  recent  review  of  experimental  data  by  Thomas  (1980)  and  a  summary 
figure  given  by  Mitra  (1979) .  The  altitude  variation  of  the  reaction 
rates  used  in  the  Mitra-Rowe  scheme  based  on  this  neutral  model  are 
presented  in  Figure  4.6. 


4.3  Steady  State  Results 

The  system  of  continuity  equations  (4.5-7)  corresponding  to  the 
chemical  model  of  Figure  4.1  have  been  solved  under  steady  state 


Figure  4.6  Altitude  variation  of  the  reaction  rates  employed  in  the  Mitra-Rowe  model 
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conditions  (i.e.  3n,/3t  =  0)  for  fixed  values  of  T  using  the  method 
±  e 

of  false  positions  as  in  the  work  of  Rowe  (1972).  Figure  4.7  shows 
the  steady  state  ion  and  electron  density  distributions  predicted  by 
the  Mitra-Rowe  scheme  for  the  daytime  ionosphere  under  unheated  con¬ 
ditions  (i.e.  T  =  T  ) .  Above  85  km  the  density  of  water  cluster  ions 
e  n 

is  small  and  equation  (4.2)  is  satisfied  while  below  85  km,  Y+  becomes 
the  dominant  positive  ion.  In  the  lower  D-region,  negative  ions  become 
important  and  one  has  [e]  «  [Y+]  =  [x']. 

Figure  4.8  shows  the  variation  of  the  steady  state  electron  den¬ 
sity  with  electron  temperature  for  the  daytime  ionosphere.  Both  n^  and 
have  been  normalized  to  their  respective  unheated  values.  Two  dif¬ 
ferent  models  of  the  electron  temperature  dependence  of  a^+  have  been 
used:  b  =  0.2  (solid  curves)  and  b  =  0.5  (dashed  curves) .  All  other 
Tg  dependent  rates  are  as  given  previously. 

At  50  km  the  change  in  n  with  T  is  controlled  by  B  with  n 

Q  0  0 

reaching  a  minimum  at  T  /T  =2.5  which  corresponds  to  the  maximum 

e  n 

value  of  B-  In  the  60  to  70  km  region,  if  one  assumes  a  slow  variation 

of  ay+  with  Te  (b  =  0.2)  then  ng  is  again  largely  controlled  by  B  • 

However,  if  a  +  varies  more  rapidly  with  T  (b  =  0.5) ,  then  the  decrease 
Y  e 

in  n  due  to  increasing  B  is  offset  by  the  increase  in  n  with  decreasing 
e  e 

otY+,  so  that  n&  begins  to  increase  before  the  maximum  value  of  S  is 

reached.  At  80  km  a  +  governs  the  variation  of  n  with  T  ,  and  n 

Y  e  e  e 

increases  more  rapidly  with  as  b  is  increased.  Above  80  km,  the 

identity  of  the  dominant  positive  ion  type  switches  from  Y+  to  NO+, 

hence  the  n  variation  with  T  is  controlled  by 

e  e  NO 


which  is  assumed 


Figure  4.7  Steady  state  ion  and  electron  density  distributions  for  unheated  conditions 
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to  vary  as  T  .  Variations  in  a  ,  effect  n  primarily  in  the 

e  02+ 

E-region  during  normal  daytime  conditions,  but  can  have  a  signifi¬ 
cant  effect  on  n^  at  D-region  heights  during  disturbed  conditions 

(Mitra,  1975b) .  Heating  induced  changed  in  a  +  have  no  signifi- 

4 

cant  effect  on  since  the  density  of  0^+  is  relatively  small. 

4 . 4  Time  Dependent  Solutions 

The  continuity  equations  for  the  six  ion  model  have  been  solved 
silutaneously  along  with  the  equations  for  electron  energy  and  power 
flow  using  basic  numerical  integration  techniques  (e.  g.  Ralston  and 
Wift,  1960).  Figures  4.9  through  4.13  illustrate  the  resultant  time 

3 

variation  of  key  ionospheric  parameters  due  to  a  period  of  10  seconds 

of  CW  heating  and  during  the  subsequent  return  of  the  plasma  to  unheated 

conditions  following  termination  of  heating.  The  results  presented 

here  are  for  heating  at  5  MHz,  X-mode  with  an  effective  radiated  power 

3  u 

of  100  MW.  The  electron  cooling  rates,  ^  e,j,  and  electron  neutral 
collision  frequency  used  in  these  calculations  were  taken  from  Banks 
and  Kockarts  (1973).  The  calculations  include  rotational,  vibrational 
and  fine  structure  cooling  as  well  as  elastic  processes.  The  ambient 
ionosphere  used  in  these  calculations  is  that  given  in  Figure  4.7.  The 
value  b  =  0.5  has  been  assumed  for  the  electron  temperature  dependence 
of  &y+ •  The  parameters  presented  in  Figures  4.9  through  4.12  have  been 
normalized  to  their  respective  ambient  (unheated)  values  which  are 
denoted  by  a  zero  subscript. 

Prior  to  the  initiation  of  CW  heating,  the  ionosphere  is  assumed 
to  be  in  a  steady  state  with  Te  =  Tn*  Once  the  heating  transmitter  is 


Figure  4.9  Electron  temperature  variation  during  and  following  a  period  of  CW  heating. 


Figure  4.10  Negative  ion  density  variation  during  and  following  a  period  of  CW  heating. 
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Figure  4.11  Electron  density  variation  during  and  following  a  period  of  CW  heating. 


Figure  4.13  Variation  of  the  negative  ion  to  electron  ratio  during  and  following 
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turned  on,  the  D-region  electron  temperature  rapidly  increases 
(Figure  4.9) ,  reaching  a  quasi-steady  state  in  about  1  ms.  Maximum 
heating  occurs  near  70  km.  In  response  to  this  temperature  enhance¬ 
ment,  the  electron  attachment  rate  increases  from  B (T  )  to  B (T  ) 

n  e 

while  the  individual  recombination  coefficients  decrease  from  a.+(T  ) 

l  n 

to  a.  +  (T  ) .  This,  in  turn,  leads  to  a  variation  in  the  ionic  constit- 
l  e 

uents . 

The  first  ion  density  change  that  takes  place  is  a  rapid  increase 
in  [02~]  due  to  the  increase  in  8  as  shown  in  Figure  4.10.  At  70  km, 
[02“]  peaks  after  1  s  of  heating,  then  decreases  as  0  “  is  converted 
to  X  .  The  net  build-up  in  negative  ions  leads  to  a  corresponding 
decrease  in  electron  density,  [e] ,  throughout  the  lower  D-region  as 
shown  in  Figure  4.11.  In  the  upper  D-region,  the  negative  ion  to  elec 
tron  ratio,  X,  is  very  small  (X<<1) ;  hence  the  electron  density  is 
unaffected  by  changes  in  the  negative  ion  chemistry.  During  the 
initial  phase  of  heating  (i.e.  t  <  10  s) ,  the  total  positive  ion  den¬ 
sity  [pj  -  [Y+],  remains  constant  (Figure  4.12).  As  [p]  gradually 
begins  to  increase  in  response  to  the  decrease  in  ay+»  [O^-]  ,  [ X~]  and 
[e]  also  increase  in  proportion  to  [p] .  As  a  consequence,  after  an 
initial  rapid  increase,  X  remains  essentially  constant,  as  shown  in 
Figure  4.13. 

Following  termination  of  CW  heating,  T  rapidly  decays  to  its 

e 

ambient  level  (Figure  4.9)  so  that  8  and  the  a^+  take  on  their  respect 
ive  ambient  values.  Thus,  [02~],  [X-]  and  [p]  decay  back  to  their 
ambient  levels,  each  on  a  progressively  longer  time  scale.  A  signifi¬ 
cant  feature  of  this  decay  process  is  that  as  the  electron  density  in 
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the  lower  D-region  recovers  from  its  depressed  level  caused  by 
heating,  it  can  overshoot  its  ambient  level  as  shown  in  Figure  4.11. 
The  reason  for  this  overshoot  phenomenon  will  become  apparent  from 
the  following  discussion. 


4 . 5  Analysis 

A  physically  meaningful  interpretation  of  the  results  presented 
above  may  be  obtained  by  considering  the  continuity  equations  for 
the  chemistry  scheme  of  Figure  4.1.  For  02~  one  has 


gf  [02~]  =  6  [e]  -  [02“]  (X  +  Y1  +  Ip])  (4.26) 

=  0  <  [p  ]  -  [X~  ])  -  [02“]  (a  +  8) 

where  a  =  X  +  y^  +  Qtj  [p].  After  the  initiation  of  CW  heating,  the 
electron  temperature  in  the  lower  D-region  quickly  attains  some  quasi¬ 
steady  state  level  and  then  remains  essentially  constant  (Figure  4.9), 
hence  8  =  8(T&)  is  fixed.  Moreover,  [x  ]  and  [p]  are  slowly  varying 
compared  to  [c>2  ],  hence  one  can  easily  solve  equation  (4.26)  and  obtain 

t°2_]  =  [°2_]  e’tAi  +  rfr  ([p]  -  [x"])(1  ■  e_tAl)  (4-27) 

o 


where 


T1  =  a  +  8' 


Similarly,  the  continuity  equation  for  [X  ]  may  be  written  as 


af  [x~]  =  X  [o2“]  -  [x~]  (y2  +  ax  [p]) 


=  C  [p]  -  [X-]  (b  +  5) 


(4.28) 
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where  b  =  Y2  +  a  [  p  ] 

X  0 

and  ?  =  rrr  • 

Again,  [p]  is  slowly  varying  compared  to  [  X  ]  and  one  may  readily  show 
that 


[X-]  =  [x']oe't/T2  +  (1  ‘  e~t/l2) 


where 


T2  =  b  +  5 


(4.29) 


The  time  constants  T  and  are  plotted  in  Figure  4.14  along 
with  the  heating  time  constant,  T  ,  and  the  time  constant  for  positive 
ion  changes,  T  ,  which  will  be  defined  presently.  It  is  apparent  from 
this  figure  that  the  time  scales  for  ion  chemistry  modifications  are 
much  larger  than  that  for  electron  temperature  changes;  hence,  one  can 
correctly  set  Sn^/Bt  =  0  in  the  energy  balance  equation  when  consid¬ 
ering  transient  electron  temperature  variations  associated  with  turning 
the  heating  transmitter  off  and  on. 

Combining  equations  (4.27)  and  (4.29f)  with  equation  (4.7),  one 
finds  that  for  t>T  >T , 

-  (1  +  S/a)  (1  +  5/b)  (4.30) 

Thus,  [p]/[e]  is  determined  by  the  ratios  0/a  and  £/b.  Typical  values 

of  a,  b,  0,  E,  and  other  parameters  based  upon  the  semi-empirical  study 

of  Rowe  et  al.  (1974)  are  given  in  Table  4.2.  The  values  are  for 

ambient  conditions.  During  heating  0  and  £  are  increased  by  some 

unknown  factor.  The  quantities  a  and  b  are  in  general  functions  of  [p] 

However,  during  the  initial  stage  of  heating  (i.e.  t  <_  T  ) ,  [  p]  =  [p], 

z  o 


region  time  constants  as  a  function  of  height 
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hence,  a  and  b  are  constant.  During  the  long  term  variation  of  [  p  ] 
(i.e.  t  >  x  ) ,  a  and  b  will  remain  constant  provided  [p]  is  small 
compared  to  and  X  +  y  .  If  one  takes  a  [p]  <  10  s  for  normal 
daytime  conditions  in  the  lower  D-region  and  X,  y  ,y2  as  given  in 
Table  4.2,  then  [p  ]<<X  +  y^  and  y^.  During  disturbed  ionospheric 
conditions  such  as  solar  flares  and  PCA's,  this  condition  may  not  be 
satisfied  (e.g.  Mitra,  1975b) .  With  a  and  b  constant  and  neglecting 
any  long  term  variation  in  6  and  one  has  by  equation  (4.30) 

[p]/[e]  =  constant  for  t  >  T  .  It  follows  that  A  =  [p]/[e]  -  1  is 
also  constant  over  this  time  scale  in  agreement  with  Figure  4.13. 

Applying  this  result  to  the  continuity  equation  for  [p]  one  has 


where 


and 


a 


V 


=  2  '  (aeff  +  “i  X) 

=  Q  -  a  [p]2 

2  =  co2+  +  2NO+ 

(C1eff  ~  <V/{(1  +  <1  +  5/b)} 


Upon  solving  equation  (4.31),  one  finds  that 


(4.31) 


[P] 


(1  -  A)  /  (1  +  A) 


(4.32) 


where  A 

and 


-  [p]0)/(«^7a  +  [p]Q)  >  e~t/TR 
1 


XR 


2  /p 


Thus,  for  a  fixed  Q,  the  magnitude  and  time  scale  for  changes  in  [p ] 
induced  by  CW  heating  are  controlled  by  a. 


The  analytic  formulation  given  above  is  also  applicable  during 
the  return  of  the  plasma  to  ambient  conditions  following  termination 
of  heating  provided  one  replaces  [ 0 ]  ,  [X  ]  ,  and  [p]Q  by  their 
respective  values  at  the  time  of  heating  termination  and  one  sets 
8  and  the  a^+  to  their  respective  ambient  values.  Thus,  the  ambient 
steady  state  positive  ion  density  is  given  by 

[p]  =  /q/5  (T  )  (4.33) 

o  n 

and  one  finds  from  (4.32)  that  the  maximum  change  in  p  due  to  CW 
heating  is 

max 

i 

and  is  attained  for  t»T  .  Referring  to  Figure  4.14,  one  can  see  that  t» 

R 

t  is  much  larger  in  the  lower  D-region  than  in  the  upper  D-region, 

K 

hence,  [p]  responds  much  more  slowly  to  heating  at  the  lower  heights. 

For  example,  in  Figure  4.12,  [p]  is  still  increasing  after  10^  s  of 
heating  at  70  km,  is  just  beginning  to  increase  at  60  km,  but  has 
reached  its  maximum  level  at  80  km. 

As  noted  earlier,  the  possibility  exists  that  Q  may  depend  on  T 

e 

and/or  E.  For  this  case,  one  can  write 


/Q(T  ,  E)a  (T  ) 
e _ n 

Q(T  )  a  (T  ) 
n  e 


(4.35) 


Thus,  the  positive  ion  concentration  can  be  enhanced  during  HF  heating 
by  decreased  recombination  or  increased  production. 


The  predicted  electron  density  variation  during  and  following 
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CW  heating  (Figure  4.11)  may  be  described  by  the  simple  formula 


«  f<T  )  W- 
e  [pi 


where  for  t>3r2. 


a  +  3  (T  )  b  +  5  (T  )  T  (T  )  T  (T  ) 
f  (T  )  =  _ 2—  - 2—  =_i - SL_ 

1  e  a  +  3  (T  )  b  +  £  (T  )  T  (T  )  t  <T  ) 
e  e  1  n  2  n 


(4.36) 


(4.37) 


During  the  initial  phase  of  heating  (t«T  )  one  has  [p]  =  [p]  ; 

hence,  the  electron  density  decreases  to  a  level  f (T^)  [e]^  on  a 

time  scale  As  [p]  gradually  begins  to  increase,  [e]  begins  to 

recover  from  its  depressed  level  on  the  time  scale  tr.  When  heating 

is  terminated,  3  takes  on  its  ambient  value,  which  results  in  a  prompt 

decrease  in  the  negative  ion  density  on  time  scale  T  2 .  On  this  time 

scale,  the  positive  ion  density  remains  fixed  at  its  heated  level, 

[p]  .  Thus  [e]  promptly  rises  to  a  level  [e]/[e]  =  [p]„/[p]  , 

ri  Otto 

resulting  in  the  overshoot  phenomena  illsutrated  in  Figure  4.11. 


Finally  on  time  scale  T  ,  [e  ]  decays  with  [p]  to  its  ambient  level. 

Two  special  forms  of  equation  (4.36)  are  worth  noting.  First, 

if  b<<£ (T  ),  one  has  f (T  )  =  3 (T  )/B(T  ),  so  that  the  prompt  decrease 
n  e  n  e 

in  e  during  heating  is  inversely  proportional  to  3.  This  condition 

may  be  satisfied  in  the  lower  D-region  if  Y2  small  (i.e.  <  0.1s 

at  60  km)  .  If  y2  is  large  (i.e.  >  4  s  '*')  the  change  in  3  due  to 

heating  may  be  obscured  so  that  f(Tg)  ^  1.  Secondly,  if  a>>B(Te>  and 

b»£  (T  )  ,  then  ot  =  o  --  and  f  (T  )  =  1  so  that 
e  eff  e 


-1 


M  ,M  .  [  y,  <V 
^  ^  II  (V  ’ 


t»T 


(4.38) 
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as  in  the  work  of  Holway  and  Meltz  (1973)  for  the  upper  D-region. 

From  the  above  analysis  it  is  clear  that  experimental  observa¬ 
tion  of  electron  density  changes  during  HF  heating  could  give  one 
valuable  information  about  the  size  of  key  ion  chemistry  parameters. 
From  long  time  scale  observations  of  [e]  during  and  following  HF 

heating,  one  can  determine  [p]max/[p]  ,  X  (T  )  and  T  (T  )  which  would 

o  R  6  R  n 

give  one  an  estimate  of  a(T  ),  a (T  )  and  Q.  Observations  of  the  prompt 

n  e 

changes  in  fe]  would  allow  one  to  determine  f  (T^)  ,  t  and  t  .  Determ¬ 
ination  of  f (T  )  could  give  one  an  indication  of  the  relative  size  of 
e 


CHAPTER  V 


HF  HEATING  INDUCED  ATTENUATION  OF  DIAGNOSTIC  WAVES 

5.1  Vertical  Incidence  Pulse  Absorption 

Measurements  of  radio  wave  absorption  by  the  vertical  incidence 
pulse  sounding  technique  (A1  method)  have  been  used  routinely  by 
numerous  experimenters  for  studies  of  the  short  and  long  term  vari¬ 
ability  of  the  ionopshere  {e.g.  Gnanalingam  and  Kane,  1973;  Ganguly, 
1974) .  The  A1  absorption  technique  consists  of  measuring  the  ampli¬ 
tudes  of  vertically  reflected,  pulsed  radio  waves  (diagnostic  waves) 
originating  from  a  ground  based  HF  sounder.  The  absorption  experienced 
by  a  diagnostic  wave  as  it  propagates  over  the  two-way  path  from  ground 
level  to  the  wave  reflection  height,  h^_,  and  back  again  is  given  in 
decibels  by  (Rishbeth  and  Garriott,  1969) 

2ui  hr 

L  =  8.686  —  f  Y  (h')  dh"  (5.1) 

c  J  aD 

h 

o 

where  u)^  is  the  angular  frequency  of  the  diagnostic  wave  and  xD  is  the 
corresponding  absorption  index.  Since  xD  depends  directly  on  the 
electron  density  and  electron-neutral  collision  frequency,  variations 
in  L  are  interpreted  in  terms  of  a  change  in  these  parameters. 

The  A1  technique  can  also  be  employed  as  a  straightforward  means 
of  monitoring  plasma  modification  resulting  from  high  power  radio  wave 
heating  of  the  ionosphere.  In  this  case,  one  measures  the  change 
in  amplitude  of  the  reflected  diagnostic  pulses  as  the  heating  trans¬ 
mitter  is  turned  on  and  off.  The  ratio  of  the  amplitude  of  the 
received  diagnostic  wave  during  high  power  heating  to  that  which  would 
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have  been  received  with  no  heating  is  given  by 

hr 

AL  =  8.686  f  [x„(h;  T  )  -  x^h;  T  )]  dh'  (5.2) 

c  D  e  D  n 

h 

Measurements  of  this  type  were  performed  at  the  Platteville  HF  heating 
facility  during  the  early  seventies  and  were  summarized  by  Utlaut 
and  Violette  (1974).  The  A1  measurements  were  made  on  a  2.667  MHz 
diagnostic  wave  and  showed  about  a  6  dB  increase  in  absorption  as  the 
plasma  was  heated  at  frequencies  between  3  and  6  MHz  with  an  effective 
radiated  power  of  about  100  MW.  The  diagnostic  wave  absorption  was 
observed  to  increase  promptly  (<  40  ms)  after  turn  on  of  the  heating 
transmitter,  but  the  time  required  for  the  absorption  to  return  to 
its  unheated  level  following  a  long  period  (10  min.)  of  CW  heating 
was  very  slow  (~  10  min.).  The  prompt  increase  in  absorption  is 
attributed  to  the  change  in  the  electron-neutral  collision  frequency 
as  the  electrons  are  quickly  heated  by  the  high  power  wave  field, 
while  the  long  term  recovery  is  most  likely  due  to  an  electron  density 
modification  arising  from  the  electron  temperature  dependence  of 
various  ion  chemistry  reactions  as  discussed  in  Chapter  XV. 

In  this  chapter,  the  electron  temperature  and  density  models 
developed  in  the  previous  chapters  are  used  to  compute  the  short  and 
long  term  variations  in  A1  absorption  due  to  HF  heating.  These  theo¬ 
retical  results  are  compared  with  the  Platteville  data  and  more  recent 
A1  absorption  measurements  at  Arecibo  by  Tomko  et  al.  (1981) . 

It  should  be  noted  that  equations  (5.1-2)  apply  only  to  the  non- 
deviative  component  of  the  absorption  experienced  by  the  diagnostic 
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wave.  In  general,  the  wave  is  also  subjected  to  deviative  absorption 
associated  with  the  bending  of  the  wave  near  its  reflection  height. 

The  reflection  heights  for  HF  waves  are  in  the  E  and  F  layers  near 
the  level  where  the  wave  and  plasma  frequency  are  equal  (i.e.  U)e=w)  . 
The  present  analysis  focuses  on  plasma  modifications  below  the  wave 
reflection  height  and  hence,  only  nondeviative  absorption  is  con¬ 
sidered. 

5.2  Collision  Frequency  Effects 

Figure  5.1  illustrates  the  change  in  absorption  as  defined  by 
equation  (5.2)  on  a  2.4  MHz  diagnostic  wave  due  to  enhancement  of 
the  electron-neutral  collision  frequency  by  HF  heating  as  a  function 
of  the  heating  frequency.  Curves  are  shown  for  X-mode  heating  (solid) 
and  0-mode  heating  (dashed)  at  100  MW  ERP  and  for  both  polarizations 
of  the  diagnostic  wave  (as  labeled) .  The  shape  of  the  attenuation 
versus  heating  frequency  curves  for  X  and  0  mode  diagnostics  are 
similar  although  in  general,  the  X-mode  diagnostic  is  absorbed  more 
than  the  O-mode  for  a  fixed  heating  frequency.  Likewise,  X  and  0-mode 
heating  for  fixed  diagnostic  polarization  result  in  similar  attenua¬ 
tion  versus  heating  frequency  curves,  although  the  peak  absorption  is 
shifted  in  frequency  between  heating  modes  due  to  the  birefrigent 
nature  of  the  ionosphere. 

The  general  shape  of  the  attenuation  versus  heating  frequency 
curve  may  be  explained  using  the  simple  quasi-longitudinal  approxima¬ 
tion  for  the  absorption  index  of  the  Appleton-Hartree  theory  (Budden, 
1966) .  With  this  approximation,  one  has 
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Heating  induced  attenuation  at  2.4  MHz  as  a 
function  of  heating  frequency. 
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where  cu  =  ou  +  u>  cos  8  (see  (2.106)).  Figure  5.2  illustrates  the 
O  f  X  D  G 

height  dependence  of  to  and  V  for  ambient  conditions,  and  the  ef- 

e  en 

fective  frequencies  w  and  u)  for  a  2.4  MHz  diagnostic  wave,  and 

U  X 

Arecibo  geomagnetic  conditions. 


Let  v  =  v  (T  ),  v  =  v  (T  )  and  Av  =  v  -  v  .  In  the  lower 
O  en  n  H  en  e  HO 

D-region,  there  are  a  range  of  heights  hQ  <  h  <  h^  (region  I)  such 

that  ci)  <  V  .  The  change  in  absorption  due  to  heating  over  this 

O  /  X  H 

range  of  heights  is  given  by 
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independent  of  the  diagnostic  frequency  or  polarization.  In  the  i5)per 

I>- region,  there  are  a  range  of  heights  (region  II)  h^<  h  <  h^  such 

that  >  V .  and  one  has; 
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(5.5) 


The  quantity  defined  by  equation  (5.4)  represents  negative  attenua¬ 
tion  or  enhancement  of  the  diagnostic  signal  and  opposes  the  effect 
of  the  absorption  term  defined  by  equation  (5.5) . 

Figure  5.3  shows  the  change  in  absorption  due  to  heating  as  a 
function  of  height  for  heating  at  3  MHz  X-mode  and  6  MHz  0-mode  at 
100  MW  ERP  and  for  both  diagnostic  polarizations.  The  enhancement 
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Figure  5.3  Heating  induced  attenuation  at  2.4  MHz  as  a 
function  of  height. 
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region  for  X-mode  diagnostic  is  clearly  defined  and  corresponds  to 
the  region  of  maximum  electron  temperature  change  (see  Figure  3.6). 

The  O-mode  diagnostic  does  not  show  a  notable  enhancement  since  the 
a)  =  v  level  is  considerably  lower  than  the  U3V  =  v„  level  which 
thus  narrows  the  range  of  integration  in  equation  (5.4) .  For  the 
lower  end  of  the  range  of  heating  frequencies  considered  here,  Av 
is  large  in  region  I  and  equation  (5.4)  is  competitive  with  equation 
(5.5).  As  one  moves  to  higher  heating  frequencies,  Av  in  region  I 
decreases  while  increasing  in  region  II,  thus  the  enhancement  term 
diminishes  and  equation  (5.5)  is  dominant.  If  only  rotational  cooling 
is  considered,  one  can  easily  show  from  equation  (3.19)  that  for  the 
higher  heating  frequencies 

3/2  2 

Av  =  C  n  .  T  '  S/(f  +  f  cos  0)  (5.6) 

axr  e  —  G 

where  C  is  a  constant  and  n  .  is  the  number  density  of  air.  Thus, 

air 

2 

AL  decreases  approximately  as  1/f  . 

U  f  A 

Comparison  of  the  theoretical  results  presented  in  Figure  5.1 
with  the  empirical  results  reported  by  Utlaut  and  Violette  (1974) 
shows  limited  agreement.  For  heating  at  3.15  MHz  X-mode  and  6.21  MHz 
0-mode  with  an  effective  radiated  power  of  100  MW,  the  observed  at¬ 
tenuation  on  a  2.667  MHz  X-mode  diagnostic  was  about  6  dB  for  both 
heating  frequencies.  The  theoretical  value  for  an  X-mode  diagnostic 
and  a  heating  frequency  of  3.15  MHz  X-mode  is  also  about  6  dB;  how¬ 
ever,  the  theoretical  value  for  6.21  MHz,  0-mode  heating,  20  dB,  is 
much  larger  them  the  observed  value.  This  discrepancy  may  be  due  in 


part  to  the  solar  zenith  angle  variation  at  the  electron  density. 

The  theoretical  values  are  for  xs  =  30°  which  is  more  appropriately 
compared  to  the  3. 15  MHz  data  which  was  taken  at  noontime  than  the 
6.21  MHz  data  which  was  taken  late  in  the  day  (1400  -  1800  LT) .  The 
variation  of  AL  with  solar  zenith  angle  is  considered  in  the  following 
section.  The  discrepancy  between  the  experimental  and  theoretical 
values  at  6.21  MHz  may  also  be  due  to  uncertainties  in  the  electron 
energy  loss  rates.  Figure  5.4  illustrates  the  effect  of  a  factor  of 
two  increase  in  the  rotational  cooling  rate.  One  notes  that  the 
theoretical  heating  induced  absorption  for  6.21  MKz  O-mode  heating  is 
substantially  reduced  with  respect  to  the  corresponding  value  without 
the  rate  change  (Figure  5.1),  while  the  value  of  AL  for  3.15  MHz 
X-mode  heating  is  not  significantly  altered. 

5.3  Electron  Density  Effects 

In  this  section,  the  effect  of  the  electron  density  distribution 
on  the  heating  induced  change  in  diagnostic  wave  absorption  is  con¬ 
sidered.  Since  the  absorbing  properties  of  the  ionosphere  are  electron 

density  dependent,  any  change  in  n  ,  whether  natural,  such  as  the 

e 

diurnal  variation  with  solar  zenith  angle  (see  Appendix  B) ,  or  heating 
induced,  can  alter  the  diagnostic  wave  attenuation.  This  is  apparent 
from  the  simple  model  for  L  given  in  the  previous  section.  Writing 

U  /  A 

out  the  electron  density  dependence  of  AL  explicitly,  one  has  from 

U  / A 


ATTENUATION  (dB) 
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Figure  5.4  Heating  induced  attenuation  at  2.4  MHz  as  a 

function  of  heating  frequency  with  the  rotational 
cooling  rate  increased  by  a  factor  of  2. 
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o 

where  A  is  a  positive  constant  and 

n  -  ambient  (unheated)  electron  density, 
eo 

An^  -  change  in  electron  density  due  to  heating, 

with  and  being  the  heated  and  ambient  collision  frequencies  as 

defined  previously.  The  first  term  on  the  right  hand  side  of  (5.7) 

describes  the  prompt  (i.e.  t  =  T  )  changes  in  AL  due  to  the  heating 

H  0 ,  X 

induced  collision  frequency  enhancement,  while  the  second  term  describes 
the  long  term  (i.e.  t  >  T  )  variation  in  AL  due  to  electron  density 

x  O  /  A 

modifications.  The  former  term  depends  on  the  ambient  electron  den¬ 
sity  which  will  vary  diurnally  as  the  electron  production  rate,  Q  , 

e 

decreases  with  increasing  solar  zenith  angle. 

Figure  5.5  illustrates  the  change  in  absorption,  AL  ,  for  a  2.4 

o 

MHz  O-mode  diagnostic  wave  during  a  period  of  10^  s  of  heating  at 
3.175  MHz,  X-mode,  120  MW  ERP,  and  during  the  subsequent  cooling  of  the 
plasma  following  heating  termination.  Curves  are  given  for  two  dif¬ 
ferent  values  of  the  solar  zenith  angle  (Xg  =  30°  and  80°)  and  two 
different  choices  of  the  electron  temperature  dependence  of  ay+ 

(b  =  0.008  and  0.5;  see  Section  4.2).  All  values  of  ALq  were  computed 
from  the  general  equation  (5.2)  using  the  Sen-Wyller  formula  (2.106) 
for  the  absorption  index.  However,  the  major  features  of  the  time 

i 

5 


I 


Figure  5.5  Heating  induced  attenuation  at  2.4  MHz  as  a  function  of  time  during  and 
following  a  period  of  CW  heating. 
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variation  of  ALq  can  be  adequately  described  by  (5.7). 

Considering  first  the  prompt  effect  due  to  collision  frequency 

enhancement  (i.e.  t  =  10  3  s  in  Figure  5.5)  one  notes  that  ALq  for 

X  =  80°  is  about  3  dB  greater  than  that  for  x  =  30°.  At  first 
s  s 

sight,  this  result  appears  to  conflict  with  equation  (5.7)  since 

o  o 

n  (x  =  30  )  >  n  (x  =  80  )  ;  however,  one  must  not  forget  that 
eo  eo  s 

any  change  in  n  will  alter  through  the  absorption  index  for 

the  heating  wave  (see  (3.19)).  Figure  5.6  illustrates  how  varies 

with  n  (x  ) •  Near  midday,  the  electron  density  in  the  middle  of  the 
eo  s 

D-region  is  relatively  large  and  the  heating  wave  is  rapidly  atten¬ 
uated  as  it  propagates  through  this  region,  leading  to  a  peak  value 
of  near  70  km.  As  Xg  increases,  n^  decreases,  causing  less  absorp¬ 
tion  of  the  heating  wave.  Consequently,  more  wave  energy  is  available 

for  heating  the  upper  ionosphere  and  the  T  peak  shifts  upwards  in 

e 

height.  From  the  analysis  given  in  the  previous  section,  one  notes 
that  there  is  more  heating  in  the  wave  attenuation  region  (i.e.  region 
II,  w  >  v  )  for  X  =  80°  than  for  X  =  30°,  while  the  heating  in  the 
wave  enhancement  region  (i.e.  region  I,  w  <  v  )  remains  about  the  same. 
Thus,  ALq  increases  with  Xg  because  the  larger  Av  in  region  II  over¬ 
compensates  for  the  decrease  in  n 

eo 

The  long  term  variation  in  AL  shown  in  Figure  5.5  follows  the 

o 

electron  density  variations  described  in  Section  4.4.  For  Xg  =  30°, 

there  is  significant  heating  in  the  lower  D-region  where  electron 

attachment  is  important,  hence  there  is  a  sharp  decrease  in  n  as  shown 

e 

in  Figure  4.11  on  a  time  scale  of  a  few  seconds.  In  this  case.  An 

e 
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is  negative,  thus  from  equation  (5.7),  Al^  decreases.  On  longer  time 

scales,  n  in  the  lower  D-region  recovers  from  its  depressed  value  on 
e 

2 

a  time  scale  of  10  s  due  to  decreased  electron  recombination  with 

hydrated  ions,  and  n^  also  increases  in  the  upper  D-region  due  to  the 

same  mechanism.  As  a  result.  Figure  5.5  shows  a  long  term  increase 

in  ALq.  As  one  expects,  the  long  term  change  in  ALq  becomes  more 

+ 

pronounced  as  the  temperature  coefficient  for  the  Y  ions,  b,  is 

o 

increased.  For  Xs  =  80  ,  most  of  the  heating  is  in  the  upper  D-region, 

hence  the  electron  attachment  effects  are  small  and  ALq  shows  only  a 

long  term  increase  due  to  the  n  enhancement  associated  with  decreased 

e 

dissociative  recombination. 

After  the  heating  is  terminated,  AL  shows  a  prompt  decrease  due 

o 

to  the  return  of  v  to  its  ambient  value  (first  term  in  (5.7)).  On 
en 

long  time  scales,  ALq  for  Xs  =  30°  exhibits  the  effect  of  electron 
density  overshoot  (see  Figure  4.11),  while  for  Xs  =  80°,  only  a  long 
term  recovery  in  ALq  is  apparent,  due  to  recombination  effects. 

5.4  Al  Absorption  Experiment 

This  section  describes  the  results  of  Al  absorption  measurements 
made  during  the  Fall,  1980,  ionospheric  modifications  program  at  the 
new  high  power  HF  transmitting  facility  of  the  National  Astronomy 
and  Ionosphere  Center  on  the  island  of  Puerto  Rico.  The  Al  absorption 
measurements  were  made  using  a  portable  2.4  MHz  ionospheric  sounding 
system  developed  at  Penn  State . 

The  NAIC  HF  heating  facility  is  located  on  the  northeast  coast  of 
Puerto  Rico  near  the  town  of  Islote,  about  50  miles  west  of  San  Juan. 
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The  facility  has  four  200  KW  transmitters  which  drive  a  4  x  8  array 
of  non-planar  log  periodic  antennas  operating  in  the  backfire  mode 
(Trask,  1979) .  The  maximum  directive  gain  of  this  array  is  expected 
to  exceed  23  dBi  over  the  3  to  12  MHz  range,  giving  an  effective 
radiated  power  of  over  160  MW.  The  half  power  beam  width  of  the  array 
is  about  10  degrees. 

The  Penn  State  A1  sounding  system  was  located  5  km  southwest  of 
the  HF  heating  facility  at  Higuillales.  This  site  was  chosen  because 
it  is  close  enough  to  the  heating  facility  that  the  vertical  sounding 
intersects  the  main  beam  of  the  heater  in  the  D-region,  it  is  in 
mountainous  terrain,  thereby  providing  good  shielding  of  the  receiving 
equipment  from  the  heater  ground  wave,  and  it  is  the  only  available 
site  where  direct  communications  with  both  the  heating  facility  and 
the  main  NAIC  facility,  Arecibo  Observatory,  were  a>.  ailable.  A  block 
diagram  of  the  sounding  system  is  given  in  Figure  5.7.  The  sounder 
consists  of  a  2.4  MHz  transmitter  with  a  peak  power  output  of  5  KW 
which  drives  an  inverted  V  dipole  antenna.  The  transmitter  is  pulsed 
at  a  rate  of  300  pulses/s  and  the  pulse  width  is  about  330  us.  The 
receiving  system  is  a  modified  version  of  that  developed  by  Sulzer 
(1973) .  The  receiver  has  a  gain  of  about  10^  and  a  bandwidth  of  about 
60  KHz  (a  compromise  between  narrow  band  for  improved  signal  to  noise 
ratio  and  a  short  rise  time  ('v  3  ps)  for  accurate  location  of  the 
detected  echo.  The  receiver  contains  automatic  gain  control  circuity 
developed  by  Sulzer  (1979)  to  keep  the  detected  receiver  output  constant 
for  wave  interaction  measurements.  The  response  time  of  the  AGO 
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system  is  about  0.05  s.  The  AGC  system  generates  a  voltage  which 
is  proportional  to  the  logarithm  of  the  RF  input  voltage  and  thus 
serves  as  an  indicator  of  diagnostic  wave  absorption.  The  voltage 
transfer  characteristics  for  this  system  are  shown  in  Figure  5.8. 

The  AGC  voltage  was  sampled  using  a  12  bit  A/D  converter  with  a 
resolution  of  5  mV.  Data  collection  and  experiment  control  were 
provided  by  a  microcomputer  and  the  digitized  AGC  voltages  were 
stored  on  5"  floppy  disks.  The  status  of  the  heating  transmitter 
(ON  or  OFF)  was  simultaneously  sampled  and  recorded. 

Table  5.1  gives  a  quick  summary  of  the  experimental  conditions 
under  which  2.4  MHz  absorption  was  measured.  On  8/29,  all  four 
heating  transmitters  were  operational  at  75  KW  on  3.175  MHz.  Based 
on  the  model  calculations  of  Trask  (1979)  the  antenna  gain  at  this 
frequency  is  about  23  dB  (allowing  3  dB  for  losses) .  The  effective 
radiated  power  for  this  day  was  thus  estimated  to  be  about  60  MW. 

On  all  of  the  other  dates  listed  in  Table  5.1,  failures  in  one  or 
more  of  the  heating  transmitters  and  their  associated  feedlines 
restricted  the  heating  capability  of  the  Xslote  HF  facility  to  two 
transmitters  operating  at  75  KW  or  less  and  driving  half  of  the  32 
element  array.  The  effective  radiated  power  on  these  dates  was 
£  15  MW.  At  this  reduced  level  of  heating,  no  correlated  variation 
in  A1  absorption  with  heater  status  (i.e.  OFF  or  ON)  was  observed. 
Accordingly,  the  data  presented  here  are  the  60  MW  results  from  8/29. 

Figures  5.9  and  5.10  illustrate  the  observed  variation  in  A1 
absorption  at  2.4  MHz  due  to  CW  heating  at  3.175  MHz  and  60  MW  ERP. 
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(SllOA)  39V110A  D9V 


Figure  5.8  Voltage  transfer  characteristic  of  2.4  MHz  receiver. 


(9P)  13A31  1VN9IS 


Figure  5.9  2.4  MHz  absorption  at  sunset  during  a  period  of  CW  heating. 


Figure  5.10  2.4  MHz  absorption  during  a  period  of  1  minute  heater  on/off  cycles. 
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The  2.4  MHz  signal  level  at  the  input  to  the  receiver  is  plotted  in 
decibels  relative  to  10  mV  as  a  function  of  time.  Typical  strong 

pp 

daytime  echos  from  an  unheated  ionosphere  are  about  10  mV  ,  while 

the  ambient  noise  level  is  <  50  uV  .  The  data  plotted  in  these 

PP 

figures  has  been  digitally  low  pass  filtered  to  a  bandwidth  of  50  mHz 
to  minimize  noise  and  short  term  variability  of  the  data  and  thus 
give  a  relatively  smooth  plot.  The  dark  bars  at  the  bottom  of  each 
figure  indicated  periods  when  the  heating  transmitter  was  on  while 
the  white  areas  indicated  period  of  heater  off. 


Figure  5.9  presents  a  23  minute  segment  of  data  at  the  start 
of  the  observation  period  on  8/29.  The  HF  heating  transmitter  was 
on  at  the  start  of  this  period  and  the  diagnostic  signal  level  was 
relatively  strong.  Sunset  occurred  shortly  after  the  start  of  the 
experiment  and  is  readily  apparent  as  a  sharp  drop  in  signal  strength 
starting  at  about  1814  LT.  A  similar  drop  in  signal  level  is  observed 
at  sunset  during  unheated  conditions.  The  heating  transmitter  was 
turned  off  at  1820  for  a  three-minute  period  and  a  rapid  increase  in 
signal  level  was  observed.  The  average  signal  enhancement  during  the 
off  period  was  8.4  +0.5  dB. 

During  the  period  1856  to  1929  LT  the  HF  heating  transmitter  was 
turned  off  and  on  at  one-minute  intervals.  Figure  5.10  shows  the  vari¬ 
ation  in  2.4  MHz  absorption  over  a  segment  of  this  period.  The  average 
increase  in  signal  level  during  the  heater  off  periods  was  about  9  dB. 

5 . 5  Discussion 

With  the  limited  new  data  given  in  the  previous  section  and  the 
preliminary  Platteville  data  reported  by  Utlaut  and  Violette  (1974), 
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one  can  form  only  a  very  crude  picture  of  the  heating  induced  vari¬ 
ation  in  diagnostic  wave  absorption  for  comparison  with  the  theoreti¬ 
cal  predictions  given  in  Figure  5.5.  The  nighttime  values  of  Al 
measured  at  Arecibo  (9  dB)  are  larger  than  the  daytime  values  obtained 
at  Platteville  (6  dB)  despite  the  fact  that  the  effective  radiated 
power  at  Arecibo  (~  60  MW)  was  probably  less  than  that  at  Platteville 
(  -  100  MW) .  This  is  consistent  with  the  diurnal  variation  of  Al 
predicted  in  Figure  5.5. 

The  slow  (10  min,)  recovery  of  AL  following  a  period  (10  min.)  of 
CW  heating  which  was  observed  at  Platteville  is  difficult  to  explain. 
The  theory  predicts  a  prompt  decrease  in  absorption  following  heating 
termination  and  then  a  much  smaller  long  term  variation  due  to  elec¬ 
tron  density  modifications.  The  overshoot  phenomenon  associated  with 
electron  attachment  partially  compensates  for  the  prompt  absorption 
decrease  as  T^  returns  to  its  ambient  level  and  it  gives  better  agree¬ 
ment  with  the  observed  long  term  recovery  of  Al  than  one  would  obtain 
if  only  recombination  effects  were  included  in  the  model. 

Finally  the  sunset  Arecibo  data  do  not  show  any  significant  long 
term  change  in  absorption  over  the  two-minute  on  and  off  periods, 
which  is  consistent  with  the  predictions  given  in  Figure  5.5. 

Evaluation  of  the  theoretical  chemistry  model  developed  in  this 
work  is  continuing  using  new  absorption  measurements  recently  taken 
at  Platteville  (Rush,  1981)  and  VLF/LF  propagation  data  (Chilton, 

1981) .  Additional  measurements  at  Al  absorption  will  probably  be 
needed  to  establish  a  broad  empirical  base  for  evaluation  and 
improvement  of  the  model. 


CHAPTER  VI 


NONLINEAR  DEMODULATION 


6.1  Basic  Theory  and  Experimental  Review 

This  chapter  shall  discuss  the  generation  of  ELF  and  VLF  radi¬ 
ations  by  HF  modification  of  the  ionospheric  current  system.  This 
phenomena  involves  the  irradiation  of  the  ionospheric  plasma  by  a 
strong  HF  wave  which  is  modulated  at  some  frequency,  The  periodic 

plasma  heating  that  results  from  this  wave  stimulates  the  emission  of 
radiation  of  frequency  Q  by  the  ionopshere;  hence  this  effect  is 
commonly  referred  to  as  nonlinear  demodulation  or  ionospheric  detec¬ 
tion. 


The  basic  physics  of  the  nonlinear  demodulation  effect  is  illus¬ 
trated  in  Figures  6.1  and  6.2.  The  electrons  within  the  area  of 
plasma  irradiated  by  the  HF  beam  experience  periodic  heating  at  the 
HF  modulation  frequency.  Since  the  ionospheric  conductivity  is  electron 
temperature  dependent,  the  conductivity  undergoes  a  similar  periodic 
variation.  Natural  ionospheric  currents  (e.g.  dynamo,  electrojet, 
etc.)  which  pass  through  the  heated  region  are  modulated  by  the  con¬ 
ductivity.  Under  ambient  conditions  the  density  of  the  ionospheric 
current  is  given  by 


J  =  0  E 
O  — O  G 


(6.1) 


where  E^  is  the  strength  of  the  geoelectric  field  which  drives  the 


currents  and  0  =  a(T  )  is  the  ambient  conductivity  tensor.  Within 

— o  —  n 


the  heated  region,  the  current  density  is 


Geometry  of  the  nonlinear  demodulation  experiment. 
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0  (T  ) 
—  e 


=  J 


Ao(T  ) 

e 


(6.2) 


where 

Aa(T  )  =  a(T  )  -  a  . 

—  e  —  e  ~o 

Thus,  HF  modulation  of  the  ionospheric  conductivity  gives  rise  to  a 
periodic  current  at  the  modulation  frequency,  fi,  of  density 


Jn  =  Ao(T  )  £  (6.3) 

This  current  is  a  source  of  ^-frequency  radiation,  a  type  of  antenna 
without  wires. 

Nonlinear  demodulation  phenomena  of  the  type  described  above 
were  observed  during  the  Soviet  heating  experiments  at  Gorkii 
(Getmantsev  et  al . ,  1974).  During  these  experiments  the  polar  currents 
were  modulated  at  frequencies  from  1  kHz  to  7  kHz  using  a  5.75  MHz 
carrier  with  an  effective  radiated  power  of  15  MW.  Signals  were  reg¬ 
ularly  observed  during  the  daytime  at  2.5  kHz  and  4  kHz  but  not  at 

_2 

1  kHz  or  7  kHz.  The  field  strength  at  2.5  kHz  varied  from  2  x  10 
-2 

to  25  x  10  UV/m.  A  distinct  diurnal  variation  of  the  signal  was 
observed  with  the  intensity  reaching  a  maximum  at  noontime  and  being 
generally  unobservable  at  nighttime.  Subsequent  measurements  by 
Budilin  et  al.  (1977)  have  shown  that  the  average  height  of  the  region 
of  generation  of  the  nonlinear  current,  J^,  is  about  75  km. 

The  above  results  were  for  magnetically  quiet  conditions. 
Getmantsev  et  al.  (1974)  also  report  that  measurements  performed  under 
magnetically  disturbed  conditions  show  that  the  VLF  field  is  highly 
variable  during  these  conditions  and  that  the  signal  strength  sometimes 
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greatly  exceeds  the  intensity  of  signals  detected  during  magnetically 
quiet  conditions.  Nighttime  measurements  by  Kapustin  et  al.  (1977) 
during  geomagnetic  disturbances  show  that  the  observed  VLF  signal 
strength  is  well  correlated  with  the  strength  and  position  of  auroral 
currents  as  detected  by  groundbased  magnetometers.  The  VLF  signals 
were  absent  during  magnetically  quiet  nighttime  conditions. 

Recently,  Turunen  et  al.  (1980)  and  Stubbe  (1980)  have  reported 
the  observation  of  VLF/ELF  signals  due  to  nonlinear  demodulation  of 
MF  and  HF  waves  during  geomagnetic  disturbances  in  Northern  Scandinavia. 
Stubbe  (1980)  observes  a  maximum  signal  strength  at  2.5  kHz. 


6 . 2  HF  Modulation  of  Plasma  Conductivities 

From  the  results  of  Chapter  III  one  can  write  the  electron  energy 
balance  equation  for  plasma  heating  by  a  modulated  HF  wave  approximately 
as  (see  (3.22)) 


_3 

3t 
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P  (t) 
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-  T  ) 
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(6.4) 


where  the  heating  production  term  is  a  periodic  function  of  the  form 


P  (t) 
e 


l 

k=-*» 


ek 


ejk&t 


(6.5) 


The  Fourier  coefficients  of  the  production  term,  P  ,  are  determined  by 

0Jv 

the  shape  of  the  modulation  envelope  on  the  heating  wave .  In  general , 

and  the  P^  are  functions  of  T^ ;  however,  in  the  following  simplified 

analysis,  they  are  assumed  to  be  constant. 

Solving  equation  (6.4)  for  T  ,  one  has 

e 

00 

T  =  T  +  Y  AT  ,  ejknt  (6.6) 

e  n  , L  ek 
k=-°° 
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where 

P  ,  T 

At  =  — SiS—  H _ 

ek  (1  +  jkftx  ) 

H 

Thus,  the  electron  temperature  exhibits  a  periodic  variation  at  ft 

and  its  harmonics.  For  kft<<l/T„,  I  AT  ,  I  attains  its  maximum  value 

H  ek' 

|p  |x  ,  while  for  kft»l/x  ,  one  has  j  AT  |  =  jp  |/kft.  Referring 

6K  il  rl  6K  SJC 

to  Figure  4.14,  one  finds  that  the  ionospheric  time  constants  restrict 
the  values  of  ft  which  can  effectively  modulate  the  electron  tempera¬ 
ture  to  the  ELF/VLF  range. 

The  periodic  temperature  variations  can  alter  the  ionospheric 

conductivity  in  two  ways.  On  time  scales  of  the  order  of  x  ,  the 

H 

collision  frequency,  v  (T  )  is  modified,  while  on  longer  time  scales 

en  e 

(i.e.  t>T^) ,  the  electron  density  will  be  modified.  The  manner  in 

which  a  depends  on  V  and  ng  is  best  illustrated  by  considering  a 

special  form  of  equation  (2.61).  If  one  assumes  that  U(v  )  can  be 

e 

replaced  by  U (T  )  (i.e.  V  iv  )  replaced  by  v  (T  ) ) ,  then  upon  inte- 

e  en  e  en  e 
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The  above  approximate  formulae  show  that  each  of  the  elements  of  the 


conductivity  tensor  has  a  distinctly  different  dependence  on  v  , 

en 
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while  the  n  dependence  is  the  same  for  all.  Electron  density  modi- 
e 

fications  to  0_  will  only  be  important  when  fd<l/T^.  Prom  Figure  4.14 

one  finds  that  SKlOO  rad/s  is  required  for  significant  electron  density 

modifications  to  occur.  For  1/T  <  ft  <  1/T  ,  only  collision  frequency 

i.  H 

effects  are  important. 


Figure  6.3  shows  the  D-region  conductivity  distribution  for 


ambient  and  heated  conditions  when  n  =  n  (T  ) ,  ft/2ir  =  1  kHz, 

e  e  n 

u)  =  5  MHz,  X-mode  at  100  MW  ERP .  The  collision  frequency  v 

en 

varies  approximately  as  T  /T  (e.g.  Banks  and  Kockarts,  1973),  thus 

e  n 


one  has  a  (T  )  <  a  (T  )  and  a  (T  )  <  a  (T  )  but  a  (T  )  >  a  (T  )  for 
He  Hn  Le  Ln  Te  Te 


M  »v  (T  )  while  a  (T  )  <  a  (T  )  for  OJ  «V  (T  )  . 
Gene  Te  T  B  Gene 


6.3  The  Wireless  Antenna 

The  principal  ionospheric  current  system  is  the  dynamo  system, 
which  is  driven  by  the  local  induced  electric  field  resulting  from  the 
movement  of  the  ionospheric  plasma  across  the  lines  of  the  geomagnetic 
field  under  the  action  of  solar  and  lunar  tidal  forces.  The  induced 
field  is  typically  of  the  order  of  1  mV/m  (Rishbeth  and  Garriott,  1969) . 
At  high  latitudes,  a  second  current  system,  the  polar  system,  becomes 
important.  Unlike  the  dynamo  currents,  the  polar  currents  are  driven 
by  electric  fields  arising  in  the  magnetosphere  (Ratcliffe,  1972) .  Both 
the  dynamo  and  polar  currents  have  peak  densities  in  the  E-region  at 
about  110  km.  The  ionospheric  currents  can  vary  considerably  during 
geomagnetic  distrubances ,  particularly  in  polar  regions,  and  they  can 
also  be  altered  by  local  sources  such  as  thunderclouds  (Park  and 
Dejnakarintra,  1973) .  Finally,  it  should  be  noted  that  recent  rocket 
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Figure  6.3  D-region  conductivity  profiles. 
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observations  suggest  that  vertical  electric  fields  of  1  V/m  or  more 
may  exist  at  times  in  the  D-region  (Kocheev  et  al.,  1976;  Tyutin,  1976; 
Hale  and  Croskey,  1979) .  If  they  exist,  these  inten.e  fields  could  be 
major  current  sources. 

As  shown  earlier,  the  natural  ionospheric  currents  can  be  modified 
by  local  conductivity  variations  induced  by  modulated  HF  heating,  giving 
rise  to  a  periodic  current  of  density  J^.  This  periodic  current  is  a 
source  of  radiation  at  the  modulation  frequency,  ft,  and  its  harmonics: 


t) 


00 


l  J 

k=-® 


kft 


(r) 


ejkftt 


(6.10) 


The  electromagnetic  field  arising  from  the  k-th  harmonic  component  of 
J  is  specified  by  the  vector  potential  (e.g.  Jordan  and  Balmain,  1968) 


t  <r  t)  =  Uo- 
Vr'  471 
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-jKR 


Jkft(r  5  I 


dr  e 


jkftt 


(6.11) 


where  R  =  |r  -  r  | .  The  range  of  integration  in  (6.11)  extends  over 
the  entire  volume  of  plasma  which  is  heated  by  the  HF  wave.  This  radiat¬ 
ing  volume  defines  the  spatial  extent  of  the  so-called  wireless  antenna. 
If  one  assumes  that  the  spatial  extent  of  the  radiating  volume  is  small 
compared  to  a  wavelength  (i.e.  |r'|<<kft/c)  then  the  vector  potential 
for  the  wireless  antenna  is  essentially  that  of  a  Hertzian  dipole  (e.g. 
Jordan  and  Balmain,  1968)  : 

V?'  fc)  ■  sfej"  Wr']  exp  •  b)  (6-12) 


Theoretical  studies  of  the  radiation  emitted  by  the  wireless  antenna  by 
Willis  and  Davis  (1973)  and  Stubbe  and  Kopka  (1977)  have  used  this 
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approximation.  The  above  formula  also  assumes  that  the  wireless  anten¬ 
na  is  surrounded  by  free  space.  In  actuality,  of  course,  the  antenna 
is  immersed  in  a  plasma,  hence  trapping  and  ducting  effects  may  play  a 
key  role  in  determining  the  vector  potential  at  a  given  point. 

6.4  The  Arecibo  HF  Demodulation  Experiment 

On  December  2  and  3,  1980,  experiments  were  performed  at  the 
Arecibo  Observatory  using  the  Islote  HF  Heating  Facility  which  sought 
to  generate  ELF/VLF  radiation  by  modulating  the  ionospheric  current  system. 
The  ionosphere  was  heated  by  radiation  from  two  transmitters  operating 
at  75  KW  each  on  a  carrier  frequency  of  3.175  MHz,  left  circular  polar¬ 
ization.  The  effective  radiated  power  was  estimated  to  be  about  30  MW. 

The  carrier  was  pulse  modulated  with  a  50%  duty  cycle  at  frequencies  of 
from  500  Hz  to  10  kHz  and  the  resultant  ELF/VLF  radiation  was  monitored 
using  the  receiving  system  shown  in  Figure  6.4.  The  receiving  antenna 
consists  of  350  turns  of  No.  14  wire  wound  on  a  plywood  form  one  meter 
in  diameter.  It  is  similar  to  that  used  by  Dinger  et  al.  (1980)  to 
measure  1-4  kHz  ambient  electromagnetic  noise.  The  signals  induced  in 
the  loop  are  amplified,  filtered  and  then  coherently  detected  using  a 
PAR  129  A  lock-in  amplifier.  Coherence  is  maintained  by  driving  the 
frequency  synthesizers  for  the  lock-in  reference  and  HF  modulation  sig¬ 
nals  with  identical  1  MHz  rubidium  frequency  standards.  The  inphase 
and  quadrature  components  of  the  detected  radiation  were  sampled  by  a 
microcomputer  controlled  data  acquisition  system  and  were  later  converted 
into  amplitude  and  phase  data.  The  lock-in  amplifier  was  set  for  a 
30  second  integration  time,  which  corresponds  to  a  system  bandwidth  of 


Figure  6.4  ELF/VLF  receiving  system. 
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8  millihertz.  The  experiment  was  generally  conducted  by  operating  the 
HF  transmitter  at  a  fixed  pulse  rate  for  2  minutes,  and  then  the  HF  was 
turned  off  for  2  minutes.  This  cycle  was  repeated  several  times,  and 
then  the  pulse  rate  was  changed  to  a  different  frequency.  The  exact 
modulation  frequencies  were  0.5,  1,  2,  2.5,  4,  5  and  10  kHz.  Strong 
signals  which  were  well  correlated  with  the  HF  status  (OFF  or  ON)  were 
detected  at  4  and  5  kHz.  The  signals  detected  at  the  other  frequencies 
were  much  weaker  and  were  not  well  correlated  with  the  HF  status. 

Figures  6 . 5  through  6 . 9  show  some  of  the  results  which  were  obtained  at 

4  and  5  kHz.  The  status  of  the  HF  heating  transmitter  is  indicated  at 
the  bottom  of  each  graph  with  the  dark  areas  indicating  that  the  HF  was 
on  and  the  white  areas  denoting  HF  off. 

Figures  6.5  and  6.6  show  the  detected  amplitude  and  phase  at  4  and 

5  kHz  respectively  during  nighttime  conditions.  The  amplitude  of  the 
detected  radiation  is  clearly  correlated  with  the  HF  status.  Amplitude 
values  are  given  in  volts  at  the  output  of  the  lock-in  amplifier.  The 
receiving  system  was  not  calibrated  during  this  preliminary  experiment; 
hence,  the  field  strength  corresponding  to  a  given  output  voltage  is  not 
known.  The  detected  signal  is,  however,  clearly  much  greater  than  the 
background  noise  level  (i.e.  HF  off) .  The  phase  characteristics  of  the 
radiation  are  well  defined  during  HF  on  and  become  random  during  HF  off. 

Figures  6.7  through  6.9  illustrate  the  diurnal  features  of  the 
detected  radiation  at  5  kHz.  At  nighttime  (Figure  6.7),  the  detected 
signal  level  was  as  large  as  12  volts.  During  the  morning  the  detected 
signal  level  was  typically  3  to  4  volts  (Figure  6.8)  ,  while  by  mid-after- 
noon  (Figure  6.9) ,  the  signal  level  was  less  than  1  volt. 


function  of  time. 


function  of  time . 


Figure  6.8  Amplitude  and  phase  of  detected  VLF  radiation  as 
function  of  time. 
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function  of  time. 


L49 


6.5  Discussion 

The  observed  correlation  between  the  detected  signal  voltage  and 
the  HF  status  could  arise  from  several  sources  other  than  nonlinear 
demodulation  effects  in  the  ionosphere.  The  most  likely  candidate 
among  these  false  effects  is  nonlinear  demodulation  of  the  HF  ground 
wave  or  the  HF  sky  wave  in  the  receiving  system.  The  observed  diurnal 
fading  of  the  detected  signal  suggests  that  its  source  is  not  the  HF 
groundwave,  since  the  groundwave  is  independent  of  ionospheric  varia¬ 
tions.  Moreover,  no  change  in  detected  signal  strength  was  observed 
when  the  loop  was  rotated.  On  the  other  hand,  the  HF  sky  wave  would 
exhibit  diurnal  features  and  could  thus  cause  the  observed  effect. 

During  24-30  March,  1981,  tests  were  conducted  at  Arecibo  to  deter¬ 
mine  whether  the  observed  signals  were  real  ionospheric  VLF  signals  or 
false  VLF  signals  arising  from  HF  demodulation  in  the  receiving  system. 
The  December  experiment  was  repeated  and  strong  signals  were  again 
detected.  During  a  period  when  the  detected  signals  were  very  stable, 
the  bandwidth  of  one  of  the  active  filters  was  adjusted  so  that  the  modu¬ 
lation  frequency  was  out  of  the  passband.  Thus,  any  real  VLF  would  be 
rejected  while  the  level  of  modulated  HF  leakage  would  remain  essentially 
the  same.  The  detected  signal  voltage  vanished  in  the  noise  when  the 
filter  was  adjusted  off  frequency,  indicating  that  HF  leakage  through  the 
filter  was  not  being  demodulated  by  nonlinear  processes  in  the  preamps 
or  the  lock-in  amplifier  to  any  significant  extent. 

During  the  March  experiments,  strong  signals  were  detected  at  all 
frequencies  from  500  Hz  to  10  KHz .  When  the  data  were  corrected  for  the 


gain-frequency  characteristics  of  the  receiving  loop,  it  was  found  that 
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the  detected  signal  strength  decreases  with  increasing  modulation 
frequency.  This  is  exactly  what  one  would  expect  from  theoretical 
consideration  since  the  change  in  collision  frequency  (or  electron 
temperature)  decreases  with  increasing  modulation  frequency  (see  equa¬ 
tion  (6.6)).  A  similar  result  is  predicted  when  electron  density 
modifications  due  to  decreased  dissociative  recombination  are  included 
(Stubbe  and  Kopka,  1977) .  It  should  be  noted,  however,  that  if  the 
primary  source  region  for  the  ELF  is  near  75  km  as  suggested  by  Budilin 
et  al.  (1977) ,  then  electron  attachment  effects  could  decrease  the 
radiated  field  strength  over  some  range  of  ELF  modulation  frequencies. 

During  the  March  experiments,  the  receiving  system  was  calibrated 
following  the  procedure  described  by  Jeans  et  al.  (1961)  using  a  single 
turn  loop  as  an  ELF/VLF  reference  source  (see  Figure  6.4).  Preliminary 
analysis  of  the  calibration  data  suggests  that  the  strength  of  the  elec¬ 
tric  field  incident  upon  the  receiving  loop  during  modulated  HF  heating 
could  be  as  large  as  20  pv/m  at  5  KHz. 

Finally,  it  is  interesting  to  note  that  the  observed  diurnal  vari¬ 
ation  of  the  detected  VLF  radiation  at  Arecibo  is  the  opposite  of  that 
reported  by  Getmantsev  et  al.  (1974)  for  polar  conditions.  At  Arecibo 
the  strength  of  the  detected  signal  is  strongest  at  nighttime  and  de¬ 
creases  as  noon  approaches,  while  at  Gorkii,  the  signal  strength  peaked 
at  noon  and  decreased  toward  nightfall. 


CHAPTER  VII 


CLOSURE 


7. 1  Summary 

This  study  has  focused  on  plasma  modifications  and  nonlinear  elec¬ 
tromagnetic  effects  arising  from  the  interaction  of  a  powerful  radio 
wave  with  the  weakly  ionized  plasma  of  the  lower  ionosphere.  The  source 
of  the  radio  wave  is  assumed  to  be  a  ground  based  HF  transmitter  with  an 
effective  radiated  power  of  the  order  of  100  MW. 

The  principal  effect  of  the  radio  wave-plasma  interaction  is  a  local 
enhancement  of  the  electron  temperature  due  to  the  absorption  of  energy 
from  the  wave  by  the  plasma  electrons.  Theoretical  calculations  pre¬ 
sented  in  this  work  show  that  the  time  scale  for  electron  temperature 
enhancements  ranges  from  a  few  microseconds  in  the  lower  D- region  to 
several  milliseconds  at  the  base  of  the  E-layer  and  that  the  steady 
state  electron  temperature  may  be  a  factor  of  five  or  more  above  the 
ambient  plasma  temperature .  The  magnitude  of  the  electron  temperature 
enhancement  is  strongly  dependent  on  the  density  of  the  plasma  electrons 
and  the  rate  of  electron  energy  transfer  to  neutral  particles  as  well  as 
on  the  frequency  and  power  density  of  the  radio  wave. 

For  a  fixed  radiated  power  and  wave  frequency,  the  radio  wave  power 
density  available  to  heat  the  plasma  electrons  at  a  given  height  is  de¬ 
termined  by  the  amount  of  absorption  experienced  by  the  wave  as  it 
propagates  to  that  height  from  ground  level.  As  the  plasma  electrons 
are  heated  by  the  radio  wave ,  the  attenuation  experienced  by  the  wave  is 
altered  since  the  absorption  properties  of  the  plasma  are  electron 
energy  dependent.  Consequently,  the  leading  edge  of  a  radio  wave  pulse 
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from  a  powerful  transmitter  can  cause  additional  absorption  on  the 
remainder  of  the  pulse  by  its  heating  action  on  the  plasma.  This  self¬ 
absorption  effect  takes  place  primarily  in  the  middle  of  the  D- region 
and  can  severely  limit  the  power  density  available  at  higher  heights. 

It  is  shown  that  the  amount  of  self-absorption  experienced  by  a  wave 
as  it  propagates  through  the  D-region  increases  as  the  radiated  power 
level  is  increased,  and  that  a  saturation  effect  may  take  place  in  which 
any  increase  in  radiated  power  beyond  some  threshold  level  is  completely 
offset  by  increased  self-absorption  resulting  in  no  net  increase  in 
power  density  above  the  D-layer.  It  is  thus  apparent  that  one  must  take 
D-region  self-absorption  effects  into  account  when  analyzing  the  results 
of  E  and  F  region  HF  heating  experiments,  and  when  one  attempts  to 
estimate  the  level  of  ionospheric  heating  due  to  microwave  radiations  by 
frequency  scaling  HF  heating  data,  as  is  done  in  ionospheric  impact 
studies  of  the  solar  power  satellite,  as  well  as  when  one  studies  D- 
region  HF  heating. 

Since  the  reaction  rates  for  chemical  processes  in  which  electrons 
participate  are  electron  energy  dependent,  radio  wave  heating  of  the 
plasma  electrons  can  modify  the  composition  of  the  plasma.  In  this 
study,  a  simplified  D-region  ion  chemistry  scheme  has  been  combined  with 
plausible  models  for  the  electron  temperature  dependence  of  electron 
attachment  and  ion-electron  recombination  processes  in  order  to  study 
plasma  density  modifications  induced  by  continuous  high  power  radio  wave 
heating.  Enhancement  of  the  electron  temperature  in  the  lower  ionosphere 
by  radio  wave  heating  causes  increased  electron  attachment  and  decreased 
recombination.  The  increase  in  attachment  produces  a  rapid  decrease 
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in  electron  density  while  enlarging  the  negative  ion  concentration  on 
a  time  scale  of  several  seconds  or  less.  On  the  other  hand,  decreased 
recombination  leads  to  increases  in  the  positive  ion  and  electron  den¬ 
sities  on  time  scales  of  several  minutes  to  an  hour.  These  two  effects 
compete  to  determine  the  electron  density  distribution  in  the  lower 
ionosphere,  with  electron  attachment  effects  being  most  important  in  the 
lower  D-region  while  ion-electron  recombination  effects  are  dominant  in 
the  upper  D-region  and  the  E-region. 

The  radio  wave  heating  induced  electron  temperature  and  electron 
density  modifications  described  above  alter  the  absorption  properties 
and  the  conductivity  of  the  ionospheric  plasma  giving  rise  to  two  import¬ 
ant  types  of  nonlinear  electromagnetic  phenomena:  the  cross-modulation 
effect  and  the  demodulation  or  detection  effect.  The  cross  modulation 
effect  is  a  two  wave  interaction  in  which  a  strong  wave  modifies  the 
absorption  properties  of  the  plasma  by  its  heating  action  and  thereby 
alters  the  amplitude  of  another  wave  which  propagates  through  the  heated 
region.  Since  the  level  of  plasma  heating  depends  on  the  amplitude  of 
the  strong  wave  field,  modulation  on  the  strong  wave  can  be  transferred 
to  the  other  wave  via  the  plasma.  In  the  detection  effect,  the  plasma 
conductivity  is  modified  in  a  periodic  manner  by  a  strong  modulated 
wave,  giving  rise  to  periodic  currents  in  the  plasma  at  the  modulation 
frequency.  These  periodic  currents  radiate  energy  at  the  modulation 
frequency;  thus,  the  plasma  effectively  demodulates  the  strong  wave. 
Theoretical  calculations  presented  in  this  work  show  that  large  changes 
in  wave  absorption  and  plasma  conductivity  are  to  be  expected  when  the 
effective  radiated  power  of  the  HF  heating  transmitter  is  of  the  order 
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of  100  MW. 

Two  experiments  were  conducted  in  order  to  observe  the  nonlinear 
electromagnetic  effects  resulting  from  HF  heating  induced  plasma  modi¬ 
fications.  The  first  experiment  used  the  vertical  incidence  pulse 
absorption  technique  to  measure  the  change  in  absorption  on  a  diag¬ 
nostic  wave  as  the  HF  heating  transmitter  was  turned  on  and  off  (i.e. 
cross-modulation  effect) .  Data  from  this  experiment  clearly  show  a 
large  (9  dB)  increase  in  wave  absorption  at  2.4  MHz  due  to  high  power 
(60  MW  ERP)  HF  heating  of  the  ionosphere .  This  result  compares  favor¬ 
ably  with  theoretical  predictions.  In  the  second  experiment,  the  HF 
heating  wave  was  modulated  at  very  low  frequencies  and  a  coherent  detec¬ 
tion  scheme  was  employed  to  observe  any  VLF  radiation  due  to  ionospheric 
demodulation  of  the  HF  wave.  .  Strong  VLF  signals  were  detected  during 
the  experiment,  and  it  is  believed  that  this  radiation  is  due  to  HF 
heating  induced  conductivity  modulation  of  the  dynamo  current  system. 

7.2  Suggestions  for  Future  Work 

At  present,  the  available  data  on  the  effects  of  high  power  HF 
heating  on  the  lower  ionosphere  and  on  radio  wave  propagation  in  this 
medium  are  still  very  limited.  There  are  a  few  observations  of  a  number 
of  different  effects  taken  at  several  HF  facilities  operating  at  various 
levels  of  heating  under  different  geophysical  conditions.  What  is 
required,  at  this  point,  is  a  systematic  study  of  HF  heating  effects 
using  several  proven  diagnostic  techniques  (e.g.  absorption,  wave  inter¬ 
action,  incoherent  backscatter,  etc.)  with  the  aim  of  determining  how 
these  effects  depend  on  the  power  density,  frequency  and  polarization  of 
the  heating  wave  and  on  changing  geophysical  conditions.  A  study  of  this 
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type  is  essential  for  testing  the  predictions  of  HF  heating  theory  and 
for  the  proper  evaluation  and  improvement  of  ionospheric  models.  From 
a  theoretical  standpoint,  several  improvements  and  extensions  of  the 
results  of  this  study  should  be  undertaken.  First,  since  the  chemistry 
of  the  D-region  is  much  less  complex  during  naturally  disturbed  condi¬ 
tions  (e.g.  solar  flares,  eclipses  and  polar  cap  absorption  events)  than 
during  the  normal  daytime,  it  would  be  useful  for  the  planning  of  future 
HF  heating  experiments  to  model  radio  wave  heating  effects  during  these 
conditions.  Secondly,  in  view  of  the  current  interest  in  the  nonlinear 
demodulation  effect  and  the  potential  applications  of  this  phenomena  to 
the  study  of  geoelectricity  and  to  submarine  communications,  a  thorough 
study  of  the  wireless  antenna  should  be  undertaken  which  incorporates 
the  plasma  heating  and  chemistry  models  developed  in  this  study  with 
empirical  and  theoretical  models  of  the  ionospheric  currents,  conduc¬ 
tivities  and  fields.  Finally,  it  should  be  noted  that  one  of  the  weak¬ 
nesses  of  the  plasma  heating  theory  presented  in  this  work  is  that  trans¬ 
port  effects  have  been  neglected;  hence,  the  importance  of  these  effects 
should  be  ascertained. 
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APPENDIX  A 

THE  ELECTRON  ATTACHMENT  COEFFICIENT 

The  three  body  attachment  coefficient  for  the  reaction 
e  +  0^  +  O^  -*•  0^  +  O^  has  been  modeled  as 


K  (T  ,  T  )  =  K  (T  )  T 
a  e  n  ao  n  e 
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where 

K  (T  )  =  1.1617  x  103  -  3.4665  T  +  3.2825  x  lcf3  T  2 
ao  n  n  n 

a  (T  )  =  7.8193  x  102  -  3.2964  T 

In  n 

a_ (T  )  =  -1.9159  x  102  +  3.7646  T  -  4.5446  x  10~3  T  2 

2  n  n  n 

a, (T  )  =  -7.5834  x  101  +  1.2277  x  10~2  T  -  7.6427  x  10_3  T  2 

3  n  n  n 

+  1.7856  x  10-5  T  3 
n 

The  constraints  used  to  construct  this  model  are: 

1.  the  electron  temperature  variation  of  the  attachment 
coefficient  for  T^  =  77  and  300°K  given  by  Chanin  et 
al.  (1962) , 

2.  the  values  of  the  attachment  coefficient  for  thermal 
equilibrium  conditions  (Tg  =  T^)  given  by  Truby  (1972)  , 

3.  location  of  the  peak  value  of  K  at  700°K  (<u>  =  0.09  eV) , 

Si 

4.  the  neutral  temperature  dependence  of  the  0.09  eV  peak 
given  by  Spence  and  Schulz  (1972) 

The  temperature  variation  of  the  model  is  illustrated  in  Figure  4.4. 
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appendix  b 

ELECTRON  DENSITY  MODELS 

The  following  table  gives  a  set  of  ambient  electron  density 
models  for  different  values  of  the  solar  zenith  angle.  The  models 
were  computed  from  the  Mitra-Rowe  six  ion  scheme  using  the  electron 
production  rate  formulae  given  by  Rowe  (1972)  . 


NOTE:  x  =  solar  zenith  angle.  Electron  densities  in  cm 
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